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ABSTRACT 

PEAPS  (Passive  Environmental  ASW  Prediction  System)  is 
a  relatively  unsophisticated  model  which  accepts  input  source 
and  receiver  parameters  and  then  predicts  sound  propagation 
characteristics  in  an  ocean  environment,  the  corresponding 
transmission  loss,  and  the  probability  of  detection.   The 
program  was  written  for  a  programmable  desk-top  calculator 
for  immediate  deployment  and  operational  testing  aboard 
small  ASW  platforms.   The  program  is  also  available  in  a 
form  suitable  for  larger  computers. 


TABLE  OF  CONTENTS 

I.  INTRODUCTION ■ 14 

A.  BACKGROUND ■ 14 

B.  PEAPS ■ ■ 15 

1.  Description — -  15 

2.  Assumptions  and  Constraints  15 

II.  SPEED  OF  SOUND  PROFILE  (SSP)  GENERATOR  17 

A.  PURPOSE — 17 

B.  INPUTS 17 

C.  MODEL — —  19 

D.  OUTPUT — —  22 

III.  RAY  TRACE  GENERATOR — 25 

A.  PURPOSE — 25 

B.  MODES ■ — 25 

C.  INPUTS — ■-— 27 

D.  MODEL 27 

E.  OUTPUT — — ~ 33 

IV.  TRANSMISSION  LOSS  GENERATOR — 38 

A.  PATH  SELECTION — 38 

1.  Surface-Duct  Path  (SDP)  — ■ 38 

2.  Reliable  Acoustic  Path  (RAP)  — 40 

3.  Bottom-Bounce  Path  (BBP)  40 

4.  Convergence-Zone  Path  (CZP)  ■ —  40 

5.  Deep  Sound-Channel  Path  (DSCP)  40 

B.  INPUTS  ■ — 40 


C.  TRANSMISSION  LOSS  MODELS 42 

1.  Surface-Duct  Path  (SDP)  42 

2.  Reliable  Acoustic  Path  (RAP)  45 

3.  Bottom-Bounce  Path  (BBP)  45 

4.  Convergence- Zone  Path  (CZP)  49 

5.  Deep  Sound-Channel  Path  (DSCP) 50 

D.  COMBINED  TRANSMISSION  LOSS — 50 

E.  OUTPUT ■ 53 

F.  SHORTCOMINGS  OF  THE  MODELS  53 

V.  PROBABILITY  OF  DETECTION  (PD)  GENERATOR  — 58 

A.  INPUTS — 58 

B.  MODEL ■ —  58 

C.  OUTPUT — ■ ■ —  59 

VI.  CONCLUDING  REMARKS — 60 

APPENDIX  A:   PEAPS  HARDWARE  DESCRIPTION ■ 61 

APPENDIX  B:   PEAPS  WORKSHEET — -  62 

APPENDIX  C:   PEAPS  OPERATOR  NOTE  SHEETS  FOR  THE 

HP9810A  -  — ■ ■ ■ 66 

APPENDIX  D:   FLOW  CHARTS,  PROGRAM  LISTINGS,  AND 

OUTPUT  FROM  THE  HP9  810A  AND  IBM  360 76 

APPENDIX  E:   PEAPS  STORAGE  LOCATIONS  IN  THE  HP9  810A  —  —  -16  3 

LIST  OF  REFERENCES — ■ — — 168 

INITIAL  DISTRIBUTION  LIST 169 


LIST  OF  TABLES 


I.  Validation  of  Sound  Speed  Parameters  20 

II.  Comparison  of  PEAPS  and  Theoretical  Range 
Increment  Determinations — ■ 34 

III.  Comparison  of  Minimum  Duct  Thickness  41 


LIST  OF  ILLUSTRATIONS 

1.  SSP  Portion  of  Operator  Worksheet 18 

2.  SSP  Printer  Output — ■ 2  3 

3.  SSP  Plot ;- — — — ■ 24 

4.  Convergence- Zone  Path  (CZP) ■ 26 

5.  Ray  Trace  Portion  of  Operator  Worksheet 2  8 

6.  Ray  Curvature  within  a  Water  Column  ■ 29 

7.  Sample  Mode  0  Run  when  CZ  Conditions  are  Met  36 

8.  Mode  2  Run  — ■ ~  37 

9.  Sonic  Layer  Depth  - — ■ — ■ —  39 

10.  Surface  Duct  Path  (SDP)  —.———.. 43 

11.  Reliable  Acoustic  Path  (RAP) — 4  6 

12.  Multiple  Bottom  Bounce  Path  (BBP)  — 47 

13.  Convergence  Zone  Gain  Model • ■ — ■ 51 

14.  Deep  Sound  Channel  Path  (DSCP) — 52 

15.  Transmission  Loss  and  Probability  of 

Detection  Plot --• ■ —  54 

16.  Comparison  of  Absorption  Coefficient  Equation 

with  Measured  Data »~ 56 

17.  Illustration  of  Merging  Techniques  for  Speed 

of  Sound  Profiles  — — ■ — —  70 


TABLE  OF  SYMBOLS 


SYMBOL 


BMA 


BO 

Bl 

GDI 

CD2 

CD3  to  CD5 


CD6 

CS1  to  CS6 

CZMAX 

CZMIN 
CZSTOP 

CZSTRT 

DFL 


DEFINITION 

The  absorption  coefficient  (dB/yd) 
derived  from  the  input  source 
frequency. 

The  angle  of  incidence  below  which  the 
bottom  loss  becomes  constant  (deter- 
mined by  the  input  source  frequency 
and  bottom  type) . 

The  intercept  used  in  the  bottom  loss 
equation. 

The  slope  used  in  the  bottom  loss 
equation. 

Surface  depth  -  the  first  point 
entered  by  the  operator  (always  0.0) . 

The  depth  of  the  sonic  layer  (if  no 
sonic  layer  exists,  then  CD2  =  CDl) . 

The  depths  of  the  critical  points 
(those  points  where  the  gradient  of 
the  speed  of  sound  changes  sign  or 
becomes  zero)  between  CD2  and  the 
bottom. 

The  bottom  depth. 

The  sound  speeds  corresponding  to  the 
critical  depths  CDl  to  CD6. 

The  maximum  range  of  the  first  conver- 
gence zone  annulus . 

The  minimum  range  to  the  first  CZ  annulus 

The  incremented  stopping  range  for  the 
CZ  segment  of  the  program. 

The  incremented  starting  range  for  the 
CZ  segment  of  the  program. 

The  diffractive  loss  out  of  the  surface 
duct;  used  in  the  surface  duct  segment 
of  the  program. 


SYMBOL 


DEFINITION 


DR 


DRPERM 


The  horizontal  range  increment  used,,, 
while  ray  tracing  =  5 [D (i) -D (i-1) ]  ' 
where  D(i)  is  the  SSP  depth  below, 
and  D(i-l)  that  above  the  ray  depth. 

Range  increment  in  the  TL/PD  program 
=  2000  yards;  used  outside  the  CZ 
annulus . 


DRWRK 


DTH 

Dl  to  D15 
G 

IBB 
IBT 


ICZ 

ICZRC 
IDSC 

IFOM 
I  PLOT 


I  RAP 


The  working  range  increment  used  by 
the  transmission  loss  program;  =  DRPERM 
outside  the  CZ  annulus,  =  0.1  the 
annulus  width  inside  the  CZ  annulus. 

The  angle  increment  to  be  added  to 
THMIN  (input  by  the  operator  when 
MODE  =  2)  . 

The  input  depths  from  the  BT  trace. 

The  sound  speed  gradient,  in  the  water 
column  segment  in  which  the  ray  is 
being  traced. 

The  bottom-bounce  propagation-path 
indicator. 

The  input  bottom  type  indicator  (1  to 
9)  which  is  converted  to  a  factor  (BT) 
used  in  the  bottom-bounce  segment  of 
the  TLOSS  program. 

The  convergence-zone  propagation-path 
indicator . 

The  CZ  range  increment  counter. 

The  deep  sound-channel  propagation- 
path  indicator. 

The  input  figure  of  merit. 

The  plotting  indicator  for  the  ray 
trace  program  used  during  MODE  =  0 . 

0  =  the  ray  trace  is  to  be  plotted. 

1  =  the  ray  trace  is  not  to  be  plotted. 

The  indicator  for  reliable  acoustic 
path. 


10 


SYMBOL 


DEFINITION 


ISD 

ISS 

IWH 
N 

NG 

NPTS 

MODE 


R 

RD 
RADINV 

RDR 
RL2 

RMAX 


The  indicator  for  surface-duct 
propagation  path. 

The  input  system  sigma. 

The  input  wave  height  (in  feet) . 

The  number  of  bottom  bounces  being 
considered  by  the  bottom  bounce 
segment  of  the  program. 

The  number  of  critical  gradient 
sections  in  the  sound  speed  profile. 

The  number  of  input  points  from  the 
BT  trace. 

The  mode  indicator  for  the  ray  trace 
program. 

0  =  automatic  ray  trace,  compute,  and 

plot. 

1  =  automatic  ray  trace,  compute  no 

plot. 

2  =  operator  enters  the  source  emission 

angles  to  be  plotted  by  the  ray 
trace  program,  then  a  Mode  i  com- 
putation to  determine  values  needed 
for  follow-on  programs. 

The  horizontal  range  of  the  ray  that 
is  being  traced  (in  ft)  and  the  TL/PD 
range  (in  yards) . 

The  input  receiver  depth  (in  feet) . 

The  inverse  of  the  radius  of  curvature 
for  a  ray  as  it  is  traced  through  a 
vertical  section  of  the  water  column. 

A  derived  factor  used  in  the  ray  trace 
program  (RADINV  *  DR) . 

The  square  of  the  cutoff  range  for  the 
reliable  acoustic  path  (RAP)  portion 
of  the  program. 

The  derived  maximum  range  to  be  con- 
sidered in  the  transmission  loss/ 
probability  of  detection  program 
(dependent  on  IFOM) . 


11 


SYMBOL 
RPT 

RS 

RT 

SD 

SE 

SF 

SS 

SSL 

STL 

SV 

SI  to  S15 

TH 

THBB 

THBL 

THMAX 

THMIN 


DEFINITION 

The  transition  range  used  in  the  deep- 
sound-channel  segment  of  the  program. 

The  derived  skip  distance  used  in  the 
surface-duct  segment  of  the  program. 

The  derived  transition  range  used  in 
the  surface-duct  segment  of  the 
program. 

The  input  source  depth,  in  feet. 

The  derived  signal  excess 
(SE  =  STL  -  IFOM) . 

The  input  source  frequency,  in 
Hertz  (Hz) . 

Speed  of  sound. 

The  derived  sea  surface  loss  used  in 
the  surface  duct  segment  of  the  program. 

The  cumulative  transmission  loss  from 
the  existing  transmission  paths  at  a 
particular  range  from  the  source. 

The  derived  speed  of  sound  at  the 
source. 

The  derived  speed  of  sound  corresponding 
to  the  input  data  from  the  BT  trace . 

The  angle  the  ray  makes  with  the 
horizontal  as  it  is  being  traced. 

The  derived  minimum  emission  angle 
from  the  source  that  will  result  in 
a  bottom  bounce  ray . 

The  angle  used  to  determine  the  bottom 
loss  in  the  bottom  bounce  segment  of 
the  program. 

The  maximum  ray  emission  angle  from 
the  source  entered  by  the  operator 
when  MODE  =  2 . 

The  minimum  ray  emission  angle  from 
the  source  entered  by  the  operator 
when  MODE  =  2. 


12 


SYMBOL 


DEFINITION 


THP 


THPP 


The  approximate  minimum  source 
emission  angle  that  will  produce  a 
convergence  zone  ray.   Used  as  a 
starting  point  in  the  automatic  modes 
(MODE  =  0  or  1)  to  determine  CZMIN 
by  iterative  ray  tracing. 

The  approximate  THBB  used  in  the  auto- 
matic modes  as  a  starting  point  for 
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I.   INTRODUCTION 

A.   BACKGROUND 

Shore-based  prediction  systems  rely  on  complex  acoustic 
models  which  require  large  computers  and  relatively  long 
computation  times.   Their  predictions  are  based  on  histori- 
cal meteorologic  and  oceanographic  data  for  regions  in  which 
environmental  conditions  are  assumed  to  be  homogeneous. 
These  systems  use  a  single  level  of  sonar  performance  as  the 
basis  for  all  range  predictions.   The  value  of  such  systems 
lies  in  their  use  as  research  tools  and  for  pre-deployment 
planning:   they  can  be  considered  to  be  strategic  systems. 

With  the  advent  of  the  programmable  desk-top  calculator 
comes  the  possibility  that  a  simple  acoustic  prediction  model 
can  be  placed  on  board  anti-submarine  warfare  (ASW)  platforms, 
This  makes  it  possible  for  the  ASW  platofrm  to  take  advantage 
of  on-scene  bathythermographic  (BT)  data  to  obtain  current 
range  predictions  based  on  the  output  of  the  platform's  own 
sensors.   In  addition,  an  on-board  system  enables  the  ASW 
platform  to  obtain  predictions  under  all  electronic  emission 
control  (EMCON)  conditions.   The  programmable  desk-top  cal- 
culators make  possible  environmental  predictions  at  the 
tactical  level. 
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B.       PEAPS 

1 .  Description 

PEAPS  (Passive  Environmental  ASW  Prediction  System) 
was  designed  for  use  in  any  programmable  desk-top  calculator 
or  general-purpose  computer  with  sufficient  capacity  for 
the  program. 

The  system  is  composed  of  four  sub-programs:   a 
speed  of  sound  profile  (SSP)  generator,  a  ray  trace  generator, 
a  transmission  loss  (TL)  generator,  and  a  probability  of 
detection  (PD)  generator.   The  SSP  and  ray  trace  generators 
are  combined  as  an  SSP/Ray  Trace  program;  the  TL  and  PD 
generators  are  combined  as  a  TL/PD  program.   These  sub- 
programs are  described  in  detail  in  subsequent  chapters. 

2 .  Assumptions  and  Constraints 

The  system  was  written  originally  in  a  machine 
language  for  use  on  a  small  capacity,  programmable,  desk- 
top calculator  (the  Hewlett-Packard  Model  9  810A  with  expanded 
storage) .   The  hardware  used  to  develop  PEAPS  is  described  in 
Appendix  A.   Total  program  capacity  of  the  calculator  was 
2036  program  steps  with  100  data-storage  registers.   This 
calculator  was  used  because  it  was  relatively  inexpensive 
and  readily  accessible;  it  was  reasoned  that  if  a  simple 
working  model  could  be  designed  on  a  small-capacity, 
programmable  calculator,  more  sophisticated  models  would  be 
possible  on  larger  programmable  calculators. 

PEAPS  is  a  relative  unsophisticated  system  which 
uses  approximation  techniques  to  compute  transmission  losses 


under  the  assumptions  of  a  flat-bottomed  ocean  and  phase 
incoherent  interference.   Lateral  variations  in  vertical 
temperature  profiles  are  neglected.   Speed  of  sound  profiles 
are  approximated  by  straight-line  segments,  with  storage 
available  for  fourteen  segments  (fifteen  input  points)  and 
five  smoothed  gradients.   When  available  to  the  user  deep 
ocean  historical  data  may  be  incorporated  into  the  input, 
but  such  data  are  neither  stored  nor  generated  internally. 
The  specific  assumptions  and  constraints  used  in  each  sub- 
program are  described  in  the  subsequent  chapter  which  deals 
with  that  sub-program. 

General  worksheets  for  PEAPS  have  been  designed  to 
guide  the  operator  through  the  sequence  of  steps  required 
to  generate  output  (Appendix  B) . 

Operator  note  sheets  also  have  been  written  in 
outline  form  to  describe  the  use  of  the  individual  sub-programs 
of  the  system  and  to  assist  the  operator  when  problems  and 
questions  arise  (Appendix  C) . 
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II.   SPEED  OF  SOUND  PROFILE  (SSP)  GENERATOR 

A.  PURPOSE 

The  speed-of-sound  generator  produces  a  plot  of  sound 
speed  versus  depth  for  inputs  from  a  bathythermograph  (BT) 
trace.   The  data  are  then  searched  for  critical  points  which 
are  stored  for  subsequent  use  by  the  ray  trace  and  transmis- 
sion loss  (TL)  programs.   A  critical  point  is  a  point  on  the 
sound  speed  profile  where  the  gradient  of  the  sound  speed 
changes  sign  or  becomes  zero.   All  sound-speed  points  are 
stored  for  use  in  the  ray  trace  sub-program. 

B.  INPUTS 

Figure  1  shows  the  SSP  portion  of  the  operator's  work- 
sheet.  The  operator  manually  inserts  the  BT  data  from  the 
worksheet  into  the  calculator  programmed  for  SSP/Ray  Trace. 
As  indicated,  the  data  entered  on  the  worksheet  are  the 
number  of  BT  points  to  be  entered,  and  the  depth  (ft) , 
temperature  (°F)  ,  and  salinity  (ppt)  of  each  input  point; 
if  salinity  is  unknown  and  historical  salinity  data  are 
unobtainable,  a  value  of  3  5  parts  per  thousand  (ppt)  is 
entered.   If  the  bottom  depth  is  greater  than  that  of  the 
BT  trace  and  historical  data  are  not  available,  it  is  recom- 
mended that  the  operator  enter  the  last  BT  temperature  as 
the  bottom  temperature.   This  will  produce  a  profile  for 
constant  temperature  from  the  last  BT  point  to  the  bottom. 
If  tables  of  historical  salinities  and  temperatures-versus-depth 
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PEAPS  OPERATOR  WORKSHEET 

(See  PEAPS  Operator  Note  Sheets) 
BT  Information  DTG: 


Location:  N/S , E/W 


1.  Turn  the  calculator  and  plotter  on.   "LOAD"  the  SSP/Ray 
trace  program  cards.   Put  the  plotting  paper  on  the 
plotter  and  set  the  "GRAPH  LIMITS"  of  the  plotter. 

2.  Press  "STOP,  END,  CONTINUE"  (this  sequence  is  used  when 
a  new  SSP  is  to  be  generated. 

Note  1 :   If  you  enter  an  erroneous  value: 

a.  If  you  have  not  pressed  "CONTINUE",  press 
"CLR  X",  then  enter  the  correct  value. 

b.  If  you  have  pressed  "CONTINUE",  you  must 
re-enter  all  the  points  again;  press  "GO  TO, 
LABEL,  A",  then  continue  with  Step  3  below. 

3.  Press  "CONTINUE". 

4.  SSP  inputs:   ENTER  number  of  BT  points  (15  maximum) 

Press  "CONTINUE"  after  each  value  is  entered  to  advance 
the  program  0utpufc 

Input       Input      Input  Speed  of  Sound 

Point  Depth  (ft)   Temp  (°F)   Salinity  (ppt)   (ft/sec) 

1        0.0       

2 


10 


15 


FIGURE  1 
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for  the  oceanic  area  of  interest  are  available  to  the 
operator,  then  use  Section  G  of  the  Operator  Note  Sheets 
(Appendix  G)  to  merge  the  historical  and  on-scene  BT  data 
for  improved  accuracy. 

C .   MODEL 

The  basic  physical  concepts  of  sound  transmission  are 
fairly  simple,  but  computation  of  the  sound  speed  for  the 
ocean  is  rather  complex  because  of  the  varying  physical 
properties  of  sea  water.   The  computation  of  sound  speed  from 
empirical  relations  based  upon  laboratory  measurements  began 
in  the  early  1950' s  (Del  Grasso,  1952;  Wilson,  1960)  . 
Attempts  to  record  accurate  hydrographic  data  while  elec- 
tronically measuring  the  sound  speed  were  successful  in  the 
early  1960's  (Hays,  1961;  Mackenzie,  1961).   For  these  latter 
experiments,  a  computed  sound  speed  was  obtained  based  on 
three  measured  parameters:   temperature,  salinity,  and  depth. 
The  results  were  then  compared  with  the  electronically- 
recorded  sound  speeds.   These  observations  resulted  in  either 
slight  modification  of  the  coefficients  in  existing  equations 
or  the  development  of  new  polynomial  equations  (Wilson,  1960; 
Leroy,  1968;  Frye  and  Pugh,  1971;  Anderson,  1971;  Del  Grasso, 
1972).   From  the  results  obtained  by  these  experimenters,  the 
variation  of  sound  speed  versus  temperature,  salinity,  and 
depth  can  be  approximately  represented  as  indicated  in  Table  I 

Approximately  9  9%  of  all  sea  water  falls  within  the 
temperature  range  -3°C  to  +30 °C  and  the  salinity  range  33  to 
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TABLE  I 
Variations  of  sound  speed  caused  by  variations  in 
temperature,  salinity,  and  depth.   A  variation  of  one 
parameter  signifies  that  the  other  parameters  are  held 
constant:   temperature  at  0°C,  zero  depth,  salinity  at 
35.00  ppt. 


Temperature  (°C) 
0-10 
10  -  20 
20  -  30 


Sound  Speed  Change 
4.5  (meters/sec) /°C 
3.0  (meters/sec) /°C 
2.0  (meters/sec) /°C 


Salinity  (ppt) 
33  -  35 
35  -  37 


1.55  (meters/sec) /ppt 
1.25  (meters/sec) /ppt 


Depth  (meters) 
0  -  10,000 


0.017  (meters/sec) /meter 
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37  ppt  (Montgomery,  1958;  Frye  and  Pugh,  1971).   Thus, 
thermal  gradients  are  usually  much  more  significant  than 
salinity  gradients  in  sound  speed  computations;  so  that  the 
temperature-depth  relation  is  a  factor  which  must  be 
considered  in  greatest  detail  in  any  sound  propagation  model. 

Because  the  Frye  and  Pugh  equation  [1]  for  sound  speed 
is  dominated  by  temperature  terms  and  requires  relatively 
few  program  steps  compared  with  other  speed  of  sound  equations, 
it  was  used  in  the  sound  speed  profile  sub-program.   The 
equation  was  modified  to  use  the  approximation 

ASS  _  ASS  AP  2.   Q  i  ASS 
AZ     AP  AZ     *   AP 

for  the  speed  of  sound  SS  in  meters/sec,  the  pressure  P  in 
bars  (gage  pressure),  and  the  depth  Z  in  meters.   The 
resulting  equation  is 


SS  =  1449.3  +  1.5848  x  10~2Z  +  1.572  x  10~7Z2 

-  3.46  x  10~16Z4  +  4.587T  -  5.356  x  10~2T2 

+  2.604  x  10"'4T3  +  1.19(S-35)  +  9.6  x  10~2  (S  -  35)  2 
+  1.354  x  10"5T2Z  -  7.19  x  10~9TZ2 

-  1.2  x  10~2(S  -  35)T 

where  SS  =  sound  speed  in  meters/sec 
Z  =  depth  in  meters 
T  =  temperature  in  °C 
S  =  salinity  in  ppt  (parts  per  thousand) . 
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The  SSP  sub-program  converts  the  input  data  to  MKS  units 
for  use  in  the  above  equation  and  then  converts  the 
resultant  sound  speeds  to  feet/sec  for  use  in  subsequent 
routines. 

The  above  equation  is  relatively  inaccurate  in  areas 
where  salinity  exceeds  36.6  ppt :   Frye  and  Pugh  claim 
accuracy  only  to  36.6  ppt  and  are  presently  working  to  extend 
the  domain  of  their  equation. 

D.   OUTPUT 

The  calculator  prints  the  input  data  and  the  corresponding 
sound  speeds  on  tape  (Fig.  2) ,  makes  a  plot  of  sound  speed 
versus  depth  (Fig.  3) ,  and  refers  the  operator  to  Operator 
Note  2  for  ray  trace  preparation  (see  Appendix  C) . 
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:.SP 

5EE  NOTE  1 


ftgure  2.  Sample  Spsod  of  Sound 
Profil©  generator 
printer  input/output. 


128. @8* 
5  8  .  6  8  * 

35. 8  8 * 

4943. 64 

17  5. 88* 
57.88 * 
35. 0 0 * 

4935.  15 


209. 8  0* 


DEPTH 
TEMP 
SflLIN 
S  P  D 

D  E  P  T  H 

TEMP 

SRL1N 

SPI! 

DEPTH 

TEMP 

SflLIN 

SPD 

DEPTH 
TEMP 
SflLIN 
SPD 


4  8  0. 88* 
58,5  8 * 
3 5.08  * 
4  8  9  8  .  3  6 

6  5  0  .  0  0  * 
4  7. 98* 

3  5.08  * 
4885. 12 

878. 8 0 * 

4  7.38  * 
35. 08  * 

4  8  8  4  .  5  6 

980. 88  * 
4  7  .  7  8  * 
35. 88  s 
4887. 7 6 

1 875. OS  * 
46. 58* 
35. 0  8 * 

4883, 3  6 


DEPTH 
TEMP 
8  fl  L  I  N 
S  P  D 

DEPTH 

TEMP 

SflLIN 

SPD 

DEPTH 
TEMP 
SHE  IN 
SPD 

DEPTH 
TEMP 
SflLIN 
SPD 

DEPTH 

TEMP 

SflLIN 

SPD 

SEE  MUTE 


1 188. 88* 

4  5  .  8  8  * 

35. O  0  * 
4  y ? 7.91 

1 3  8  8. 8  0 * 
45. 08  * 

3  5  .  8  8  * 


1388. 00* 

.-■  A  C   l-J  v 

H  "4-  =  -J  fcl  w 

35.0  8  * 

4873. 2  6 

1580.  8  8  * 

4  3  .  7  0  * 
3 5 . 80 * 
4  8  6  9  .  4  9 

75  0  0.  8  8  £ 
43.78* 

35. 8  0  * 

4  9  6  6 . 7  8 
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III.       RAY    TRACE    GENERATOR 

A.  PURPOSE 

A  ray  trace  sub-program  included  in  PEAPS  makes  the 
following  tests  and  decisions: 

1.  inspect  the  speed-of-sound  profile  to  determine 
whether  the  characteristics  necessary  for  convergence  zone 
propagation  are  present  or  not; 

2.  if  CZ  conditions  exist,  determine  the  maximum  and 
minimum  CZ  ranges,  where  minimum  CZ  range  is  defined  as  the 
range  from  the  source  to  the  inner  limit  of  the  CZ  annulus 
and  maximum  CZ  range  is  defined  as  the  range  to  the  outer 
limit  of  the  CZ  annulus  (Fig.  4) ,  and  in  any  event; 

3.  determine  the  minimum  bottom  bounce  angle. 

B.  MODES 

The  ray  trace  sub-program  has  three  modes  of  operation 
designated  by  0,  1,  and  2. 

Mode  0  automatically  determines  and  plots  the  limiting 
rays  of  the  convergence-zone  ray-family  and  stores  the 
corresponding  ranges. 

Mode  1  is  the  same  as  Mode  0,  but  omits  the  ray  trace 
plot. 

Mode  2  gives  a  plot  of  the  operator's  input  ray-family 
followed  by  a  Mode  1  run. 
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C.  INPUTS 

The  operator  enters  the  Mode  number  and  source  depth 
from  the  ray  trace  portion  of  the  worksheet.   If  Mode  2 
is  selected,  he  must  also  enter  twice  the  maximum  and 
minimum  emission  angles  and  twice  the  desired  increment  in 
emission  angle  between  successive  rays   (Fig.  5) . 

D.  MODEL 

The  ray  trace  generator  uses  Snell's  Law  in  a  small 
angle  approximation. 

Snell's  Law  for  rays  can  be  written  as 

Ar  =  R(sinG2  -  sin01) 
and 

AZ  =  R(cos61  -  cos62) 

where  the  geometry  is  suggested  in  Fig.  6.   Elimination  of 
82  gives  the  range-depth  equation  in  terms  of  the  initial 
angle  of  depression, 


AZ1  =  Rcos91  -  R[l  -  (sinOl  +  ~)2]^  (1) 

R 


See  Section  B  of  Note  2  in  the  Operator  Note  Sheets 
(Appendix  C) . 


27 


Note  2 ;   Change  the  plotting  paper  on  the  plotter.   Enter 
the  ray  trace  program  MODE  indicator. 

0   =  Automatic  compute  and  plot 
M0DE:  ______  1  =  Automatic  compute,  no  plot 

2  =  Manual  inputs  for  plotting,  Automatic 
compute . 

If  a  ray  trace  input  error  is  made,  press  "GOTO, 
LABEL,  B,  CONTINUE",  then  input  the  MODEL  indicator 
and  ray  trace  data  again.   (This  sequence  is  used 
to  generate  a  new  trace  using  the  old  SSP  data) 

Press  "CONTINUE" 

4.  P.ay  trace  inputs: 

Press  "CONTINUE"  after  each  input  to  advance  the  program, 

Source  Depth: (ft) 

If  Mode  2  is  used: 

Minimum  Ray  Angle  x  2:  (degrees) 

Angle  Increment  x  2:    (degrees) 

Maximum  Ray  Angle  x  2:  (degrees) 

Note  3 :   If,  while  in  Mode  0   and  after  the  bottom  limiting 
ray  has  been  drawn,  a  channelled  ray  is  drawn  to 
the  plot's  range  limit;  see  Operator  Note  sheet 
section  Note  2E . 

Change  the  plotting  paper  on  the  plotter.   Press 
"STOP,  END".   LOAD  the  TLOSS/PD  program  cards. 
Press  "STOP, END".   (The  "STOP,  END"  sequence  is 
used  when  a  TLOSS/PD  input  error  is  made,  or  a  new 
TLOSS/PD  plot  is  to  be  generated  based  on  the  old 
SSP/Ray  trace  data.) 

5.  TLOSS/PD  inputs: 
Press  "CONTINUE"  after  each  input  to  advance  the  program. 

Press  "CONTINUE" 

Source  Frequency: __(Hz) 

Figure  of  Merit  (FOM)  : (dB) 

System  Sigma:         (use  10  if  unknown) 

Wave  Height:  0   1  2  3  4  5  6  7  8  9  10  11  12  (ft) 

Receiver  Depth:        (ft) 

Bottom  Type:    123456789 

FIGUPvE    5 
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•isSfi) 


/\r  =  r2  -  rl 
fcz  *  z2  -  zl 


figure  6,      Illustration  of  ray  curvature  within 
a  water  column. 
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Expanding  Azi  as  a  MacLaurin  series  about  Ar  =  0  gives 


AZ2  =  tanei(Ar)  +  1         -  (Ar)2  +  —  Sln61 =■'  (Ar)  3 

2R(cos01)  2R  (cos01)b 


+  1  -  4(sin61)2  (Ar)4  +  (2) 

8RJ(cosei) 


Substitution  of  the  trigonometric  series  representations 
of  tan  x  and  cos  x  into  the  first  two  terms  of  equation  (2) 
gives 

AZ3  =  eld  +  li_)  (Ar)  +  n_J_i^iLJ_  (Ar)2  m  (3) 

To  conserve  program  space  and  reduce  the  computation 
time,  PEAPS  uses  a  simplified  form  of  Eq.  (3) .   Use  of  the 
small  angle  approximations 

sine  =  9 

and 

e2 

COS0  »  1  -  y  , 

v/here  6  is  in  radians,  yields  the  simple  approximations 


2 

AZ  =  (61)  (Ar)  +  (A^  (4) 


and 

AG  =  Ar/R  . 
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'or  Eq.  (4)  to  accurately  represent  Eq.  (3) ,  restrictions 

o 

mst  be  imposed  on  Ar  and  91.   In  order  that  omitting  01 

ntroduces  no  more  than  10%  error,  61  must  be  less  than 
.33  radians.   To  determine  the  limits  on  Ar,  examine  the 
ase  of  "worst  possible"  conditions: 

1.  A  maximum  gradient  of  0.5  sec  and  a  minimum  sound 
>peed  of  4  70  0  ft/sec  are  assumed  so  that  the  minimum  radius 
>f  curvature  is  9400  ft. 

2.  The  maximum  allowable  angle  of  depression  of  0.35 
radians  (20°)  is  taken. 

'o  ensure  that  no  layer  is  missed  or  poorly  represented  in 
he  ray  trace ,  five  depth  increments  are  required  within  any 
.ayer  of  depth  D.   Then,  the  first-order  term  in  Eq.  (4) 
'ields  the  equation 


Ar  <   5—  =  0.57D  (5) 


is  an  approximate  bound  on  the  allowed  Ar  within  a  layer 
)f  depth  D. 

So  that  the  truncation  of  Eq.  (2)  to  second  order  in  Ar 
Introduces  no  greater  error,  the  third  term  of  the  equation 
nust  be  much  less  than  the  second, 


7sin61   ,-  (Ar)3  «  i -,  (Ar)2 

2!T(cosei)  2R(cosel) 


Ar  «  Resell!  <  2.6R 

smGl    = 
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If  Ar  is  required  to  be  at  least  one  order  of  magnitude 
smaller  than  2.6R,  then 

Ar  <_  0.26R  <    2444  ft. 

However,  since  for  small  91  the  third  term  becomes 
arbitrarily  small,  it  is  also  necessary  to  require  that  the 
fourth  term  of  Eq.  (2)  be  much  less  than  the  second, 


1  +  «(°inei)2  (Ar)4  « 1 _  (Ar)2  , 

8RJ(cos61)  2R(cos61)J 


or 


Ar  <<  (2R)  (cosOl)2  [1  +  4(sin61)2p  =  2.1R  . 


If  Ar  is  to  be  at  least  one  order  of  magnitude  less 
than  2.1R,  then 

Ar  <    0.21R  <  2012  ft. 

To  estimate  the  error  in  calculating  AZ  from  Eq.  (4) , 
subtract  Eq.  (4)  from  Eq.  (3): 


3  2 

AZ3  -  AZ  =  2±-  (Ar)  +  ~-    (Ar)2  (6) 
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The  fractional  error  is  given  by: 


AZ3  -  AZ    0.014  +  (9.56  x  10"6) (Ar) 


AZ3       0.36  +  (1.16  x  10  4) (Ar) 


Away  from  a  turning  point,  for  large  angles  where  the 
difference  defined  by  Eq.  (6)  is  large,  the  error  is 
approximately  4%  for  small  Ar  and  approximately  8%  for  large 
Ar.   If  the  error  is  required  to  be  no  greater  than  5%,  then 
Ar  must  be  no  greater  than  about  1000  feet  (near  a  turning 
point,  AZ3  -  AZ  approaches  zero  very  rapidly,  so  that  errors 
in  depth  are  negligible  for  small  angles  of  elevation  or 
depression) . 

Based  on  these  error  requirements,  the  simple  equation 


Ar  =  5(D)** 


was  used  for  Ar  and  D  in  feet.   For  the  ranges  of  D  employed 
in  PEAPS  the  restrictions  on  Ar  set  forth  above  are  not 
seriously  violated  (Table  II) . 

E.   OUTPUT 

If  convergence  zone  (CZ)  conditions  are  not  present, 
"NO  CZ"  appears  on  the  printer.   In  all  cases  "SEE  NOTE  3" 
and  "END"  appear,  indicating  the  end  of  the  SSP/P.ay  Trace 
programming  sequence,  and  referring  the  user  to  Operator 
Note  3  (Appendix  C) . 
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TABLE  II 
Comparison  of  PEAPS  and  theoretical  range  increment  values 


D 

Ar  = 

=  min(0.5  7D, 
28.5 

1000) 

Ar  =  5[D]^ 

50 

35.3 

100 

57.0 

50.0 

200 

114.0 

70.7 

400 

228.0 

100.0 

800 

456.0 

141.4 

1600 

912.0 

200.0 

3200 

1000.0 

282.8 

6400 

1000.0 

400.0 

Graphic  outputs  depend  upon  the  Mode  number  chosen. 
figure  7  is  an  illustration  of  a  Mode  0   run  when  CZ  conditions 
are  met.   Figure  8  shows  a  Mode  2  run.   These  pictorial 
outputs  may  serve  as  training  aids  in  either  a  classroom 
Dr  a  shipboard  environment  to  demonstrate  the  dependence  of 
:ay  paths  on  the  vertical  sound  speed  structure.   Mode  1 
las  no  graphic  ray  trace  output. 
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IV.   TRANSMISSION  LOSS  GENERATOR 

A.   PATH  SELECTION 

Transmission  loss  is  computed  for  paths  which  meet 
the  following  criteria: 

1.   Surface  Duct  Path  (SDP) 

a.  The  sound  speed  profile  (SSP)  contains  a  sonic 
layer.   A  sonic  layer  is  defined  as  the  layer  above  the 
first  negative  sound-speed  gradient  (Fig.  9) . 

b.  The  source  and  receiver  are  both  located  within 
the  sonic  layer. 

c.  Sound  remains  trapped  in  the  sonic  layer.   This 
trapping  fails  to  occur  at  very  low  frequencies,  where  the 
wavelength  has  become  too  large  to  "fit"  in  the  duct.   Using 
a  sound  speed  of  5170  ft/sec  and  a  sound  speed  gradient  of 
0.017  sec   in  the  sonic  layer,  we  have  from  Urick  [2] 


H  .   -  36.64  A0*667, 
mm 


for  the  minimum  duct  thickness  (H  .  )  in  feet  which  will 

man 

trap  sound  with  wavelength  A  in  feet.   Urick  states  that 
this  expression  does  not  represent  a  sharp  cutoff.   To 
reduce  programming  steps  for  the  calculator,  the  H  . 
equation  was  simplified  to 


H  .   -  32  A0'7 
mm 
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CD1 


CD2 


C(z) 


CD1 


CD2 


C(z) 


p-C(s) 


CD1  =  Surface 

CD2  =  Sonic  Layer  Depth  (SLD) 


figure  9.   Illustrations  of  Sonic  Layer  Depths  (SLD). 
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Table  III  compares  the  frequencies  and  corresponding  minimum 
duct  thickness  as  determined  by  PEAPS ,  Urick,  and  ASRAP  III 
correction  tables. 

2 .  Reliable  Acoustic  Path  (RAP) 

RAP  conditions  always  exist  except  when  the  source 
and  receiver  are  both  located  in  a  sonic  layer. 

3.  Bottom  Bounce  Path  (BBP) 
Bottom  Bounce  Path  always  exists. 

4 .  Convergence  Zone  Path  (CZP) 

CZP  is  computed  if  an  SSP  with  a  sound  speed  excess 
is  present  (Fig.  4)  . 

5.  Deep  Sound  Channel  Path  (DSCP) 

The  Deep  Sound  Channel  Path  exists  when  the  source 
and  receiver  are  both  located  in  the  deep  sound  channel  and 
the  sound  speed  profile  can  support  CZ  propagation. 

B.   INPUTS 

The  operator  loads  the  TL/PD  program  and  enters  (from 
the  TL/PD  portion  of  the  worksheet) (Fig.  5)  source  frequency 
(Hz) ,  Figure  of  Merit  (dB) ,  System  Sigma  (dB) ,  wave  height 
(ft) ,  and  MGS  bottom  type.   Although  the  figure  of  merit 
ind  system  sigma  are  used  only  in  the  PD  generator,  they 
ire  entered  at  this  time  to  simplify  operator  procedures. 

The  transmission  loss  sub-program  uses  the  following 
values  from  the  ray  trace  generator:   a  CZ  indicator,  the 
\aximum  and  minimum  CZ  ranges  where  CZ  conditions  exist, 
.he  minimum  bottom  bounce  angle,  and  the  depths  and  speeds 
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TABLE  III 
Comparison  of  minimum  duct  thickness 


Frequency 
(Hz) 

(ft) 

53 

97.5 

104 

49.7 

208 

24.9 

588 

8.8 

1080 

4.8 

3055 

1.7 

Minimum  Duct  Thickness  -  Hm^n  (ft) 
PEAPS        Urick 


789.7  776.2 

492.7  495.3 

303.7  312.4 

146.6  156.2 

95.9  104.3 

46.4  52.2 


ASRAPS  III 

750 

475 

300 

150 

100 

50 
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of  sound  at  the  critical  points  of  the  sound  speed  profile 
(sonic  layer,  deep  sound  channel  axis,  etc.). 

C.   TRANSMISSION  LOSS  MODELS 

The  program  looks  at  each  path  selected  and  determines 
transmission  losses  according  to  certain  simple  models. 
Names  used  in  this  section  are  those  used  in  the  FORTRAN 
version  of  the  REAPS  program  (Appendix  D) . 

1  •   Surface  Duct  Path  (SDP) 

Figure  10  illustrates  SDP. 


TL0__,  =  201og(R)  +  AR  +  DFL*R   for  R  <  RT  ,  R  <  RS 


TL0_D  =  101og(R)  +  AR  +  DFL*R  +  lOlog(RT) 
for  RT  <  R  <  RS 


TLSDP  =  10  lo9(R)  +  AR  +  &FL*R  +  lOlog(RT)  +  SSL*R 

for  R  >  RS  ,  R  >  RT 


TLSDP  =  201o<?(R)  +  AR  +  DFL*R  +  SSL*R   for  RS  <  R  <  RT 


where  R  =  the  horizontal  distance  from  the  source  to  the 
receiver  in  yards 

A  =  the  absorption  coefficient  [3] , 


=  [£^4  ,  -_i°Ei7]  *  ID"3   (dB/yd)  , 
1+F     410  0  -iF 


F  =  frequency  in  kilohertz  (KHz)  ,• 
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DFL  =  diffractive  leakage  out  of  the  duct  based  on 
Normal  Mode  Theory  [4  and  5],  given  by 


744f"1,66g"0,33(CD2)"3  *  10"3  (dB/yd)  ,  and 

CS3  -  CS° 
G   =  the  below- layer  gradient  =  prn  -  cn2  ' 

CD2  =  the  sonic  layer  depth  (ft) ; 

CD3  =  the  depth  below  CD2  at  which  the  next  positive 
gradient  begins  or  the  bottom  depth  (ft) ; 

CS2,  CS3  =  the  sound  speeds  (ft/sec)  which  correspond 
to  CD2  and  CD3,  respectively 

RT   =  the  transition  range  (the  range  at  which  the  duct 

is  considered  full  and  cylindrical  spreading  begins) 

,   r~,  CD2 

in  yards  [2]  =  -? /g„i-tf0-r 

-*  6arccos  (SV/CS2) 

SV  =  the  speed  of  sound  at  the  source  (ft/sec) , 

RS   =  the  skip  distance  (the  range  between  consecutive 
surface  reflections  of  the  deepest  trapped  ray) 

,   ro,    2arccos (CS1/CS2)  ,  ,, 
in  yards  [2]  = 3(R^DIN^ (yd) 

RADINV  =  the  inverse  of  the  radius  of  curvature  of  a 
ray  within  the  surface  duct  =  GD/SD. 

SD  =  the  depth  of  the  source  of  sound  (ft) . 

GD  =  the  effective  gradient  in  the  sonic  layer 

CS2  -  CS1 

CD2  -  CDl 

0  73 
0  OB  * 
SSL  =  sea  surface  loss  =  — ~ (dB/yd) 

Kb 

B    =  F*H 

H   =  wave  height  (ft) . 

The  expression  for  SSL  was  developed  from  a  power  regression 
analysis  of  data  presented  in  Reference  6,  and  is  valid  for 
values  of  B  between  0  and  30. 
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2.   Reliable  Acoustic  Path  (RAP) 
From  Fig.  11, 


TLRAP  =  201°9(K>  +  A*K  ' 


where  K  is  the  slant  range  from  the  source  below  the  layer 
to  the  surface  at  range  R,  in  yards. 

The  maximum  RAP  range  =  RL  +  RS/2 ,  where  RL  is  the 
cutoff  range,  defined  as  the  horizontal  distance  from  the 
source  to  the  point  at  which  a  ray  just  grazes  the  bottom 
of  the  sonic  layer  (or  surface  if  no  sonic  layer  exists) , 
and  RS  is  the  skip  distance  (described  previously  in  the 
surface  duct  model) . 

3 .   Bottom  Bounce  Path  (BBP) 

Four  separate  bottom-bounce  paths  are  possible  between 
the  source  and  receiver:   two  with  one  surface  reflection, 
one  with  two  surface  reflections,  and  one  with  no  surface 
reflections  (Fig.  12) .   Incoherent  superposition  of  these 
rays  at  the  receiver  results  in  a  6  dB  reduction  in  the 
transmission  loss.   For  emission  angles  greater  than  about 
15° ,  the  radius  of  curvature  of  a  ray  is  usually  large  enough 
to  assume  straight-line  propagation  of  the  rays  from  the 
source  to  the  bottom  and  back  to  the  surface.   A  ray  from  the 
source  will  strike  the  bottom  at  an  angle  of  incidence 
estimated  by 
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mITnT       -lrCD6  -  SD, 
THBL  =  tan   [ r= ] 


where 


R  (CD6  -  SD)   f  .  . 
XB    N  (2CD6  -  SD)  lYdS) 


is  the  horizontal  distance  between  the  source  and  the  range 
of  the  first  bottom  bounce,  CD6  is  the  bottom  depth  (ft) , 
and  N  is  the  number  of  bottom  reflections. 

The  distance  traveled  by  the  four  rays  which  strike 
the  bottom  once  can  all  be  approximated  by 


(yds) 


cos (THBL) 

for  reasonably  small  source  and  receiver  depths  or  large 
range.   Therefore,  an  approximate  transmission  loss  can  be 
written  as 


TLDT3  =  201og[ ,?„T.Ti  ]  +  — ^L1^-~—  +   n*BL 

BB        ^  Lcos (THBL)  J    cos (THBL) 


where  the  subtraction  of  6  takes  into  account  the  6  dB  gain 
caused  by  multiple  bottom  bounce  paths,  and  BL  is  a  bottom 
loss  term  which  was  determined  through  regression  analysis 
of  bottom  loss  data  presented  in  Reference  7.   Specific  values 
for  BL  are  determined  by  a  sorting  routine  which  uses  the 
input  frequency  and  bottom  type.   THBL  is  the  angle  at  which 
the  ray  strikes  the  bottom. 
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4 .   Convergence  Zone  Path  (CZP) 
See  Fig.  4. 


TL_„  =  201og(R)  +  A*R 

cz 


when  R  is  outside  the  CZ  annulus  and 


TL__  -    201og(R)  +  A*R  -  GCZ 


when  R  is  bounded  by  CZSTRT  and  CZSTOP.   The  CZ  annulus  is 
defined  as  the  width  of  the  zone,  CZSTOP  -  CZSTRT,  in  yards. 
GCZ  is  the  reduction  in  transmission  loss  caused  by  the 
increase  in  sound  pressure  level  (SPL)  which  results  from  ray 
focusing  within  the  ray  family  forming  the  CZ  annulus. 


PP  9 

GCZ  =  ~  exp(-RPV2GB) 


»  =¥^^H 


r*r>  0  .  36      .   ,-2. 

gb  =  — - — y~   (yd  ) 

GCZMAX 


RANGE  =  R  -  CZMIN  (yds) 

X   =  CZ  annular  width  (yds) 

GCZMAX  =  15.454  5  -  4.54  5  *  10_4X  (yds) 

The  equations  for  GCZ  and  GCZMAX  are  empirical  and  were 
derived  from  analysis  of  numerous  transmission  loss  curves 

9 

obtained  as  output  from  the  FACT  model  [5] .   The  analysis 
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established  that  the  gain  distribution  within  a  convergence 
zone  closely  resembled  a  Rayleigh  distribution  and  that  the 
gain  tended  to  a  maximum  at  approximately  18%  of  the  annular 
width.   Thus,  GCZ  was  modelled  as  a  Rayleigh  distribution 
with  the  peak  occurring  at  a  distance  from  CZSTRT  equal  to 
18%  of  the  annular  width  (Fig.  13) . 

5.   Deep  Sound  Channel  Path  (DSCP) 
See  Fig.  14. 


TLDSCP  =  101og(R)  +  A*R  +  lOlog(RTP) 


where  RTP  is  the  deep  sound  channel  transition  range, 


RTP  =  CZMIN  {D, h 


D   =  CD3  -  CD2  in  feet, 

CD3  =  depth  of  the  sound  channel  axis  (ft) , 

Z    =  SD  -  CD2  (ft) , 

CZMIN  is  in  yards. 

The  transmission  loss  term  is  valid  only  for  R  >  RTP, 
although  it  is  used  for  all  R  when  this  path  is  present. 

D.   COMBINED  TRANSMISSION  LOSS 

Since  transmission  loss  caused  by  bottom  bounce  is  always 
present,  the  total  transmission  loss  (STL)  is  initialized 
with  TL   after  each  increment  of  range.   Transmission  loss 
for  subsequent  propagation  paths  meeting  their  respective 


50 


figure   13.     Convergence  Zone  Gain  Model 


GCZMAX 


0.18(CZST0P   -   CZSTRT) 


KP 


GCZMAX    =   15.5   -    (0. 000^5 (CZSTOP   -   CZSTRT)) 

GB  =  0,36  /  GCZMAX2 

RANGE  =  r  -  CZSTRT 

RP  =  (RANGE  xVgB) 


TomrrczsTOP  ~  czstrtTI 


GCZ  =  /RP] exp 
\GB 


-RP 
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* 


existence  criteria  are  computed  for  the  same  range  and 
combined  with  STL  as  follows:   If  ©  signifies  the  combining 
of  decibels  (dB)  for  independent  (phase  incoherent)  signals, 
then 

STL  ©  TL(i)  =  STL  -  A 

where  TL(i)  is  the  transmission  loss  for  the  specific 
propagation  path  (CZ,  SDP,  DSC,  RAP ,  as  they  exist). 


A  =  101og(l + 10"0,1J) 
J  =  STL  -  TL(i) . 


E.  OUTPUT 

The  program  plots  the  combined  transmission  loss  (STL) 
in  dB  versus  range  in  nautical  miles  (NM) ,  with  a  reference 
line  for  the  input  figure  of  merit  (Fig.  15). 

F.  SHORTCOMINGS  OF  THE  MODELS 

Since  the  models  were  developed  for  use  in  a  relatively 
small  programmable  desk-top  calculator,  limitations  in  the 
available  storage  and  the  number  of  program  steps  have 
resulted  in  certain  shortcomings: 

1.  Phase  coherence  is  neglected. 

2.  Sea  surface  scattering  is  not  taken  into  account  in 
calculating  TL  for  the  convergence  zone  model. 

3.  If  the  source  is  in  the  sonic  layer  and  the  duct 

thickness  is  less  than  H  .  .  propagation  in  the  layer  is 

mm        J 

neglected. 


5  3 


t r~ — r 

+ 

+ 


T" 

+ 

+ 

+ 

+ 
+ 


N 


n 

I 

Z 
u 

hi 
ID 

Z 

c 


h — + 


+#  + 


+ 

+ 
+ 

+ 
+ 

+ 


© 


© 
s> 
m 


© 
© 


© 

CO 


© 
© 


© 
3 


© 
© 

o 


+ 


"3    d    d    d~ 


© 
© 


H 

w 

0< 

d) 

E-" 

W 
Q 

£ 

?f 

§ 

pj 
> 

En 

o 

H 

rtj 

y 

(J 

Cn 

U} 

£» 

« 

W 

o 


0) 
H 

P<EH-P 

e  w  o 

WO<  ft 


-       +    * 

W 

•H 


+ 

+ 

+ 
+ 


+ 
+ 
+ 


_| J_| J. ( 1 f__J J f 

OlfahlflblTPiN-B 

BHSHBEHHSfs 

i30d 


bsbkhe~npi 

MIOhlDOl---- 
tS<J]  55D1J. 


54 


4.  The  use  of  straight-line  propagation  will  induce 
some  error  in  the  angles  computed  for  bottom  bounce 
parameters  and  in  the  transmission  loss  computation. 

5.  Convergence  zones  are  considered  when  any  sound 
speed  excess  is  present. 

6.  The  convergence  zone  gain  term,  GCZ ,  is  found  from 
an  empirical  formula:   there  is  no  physical  justification 
for  the  use  of  the  Rayleigh   distribution. 

7.  The  model  for  the  deep  sound-channel  overestimates 
transmission  loss  for  ranges  less  than  the  (deep  sound- 
channel)  transition  range  . 

8.  For  reliable  acoustic  path,  steep-gradient  curvature 
is  not  considered;  transmission  loss  is  computed  for  straight 
line  distances  between  the  source  and  receiver. 

9.  The  absorption  loss  term  uses  Thorp's  formula  for 
all  frequencies.   As  shown  in  Fig.  16,  the  formula  is  good 
for  frequencies  greater  than  2  0  Hertz. 

10.   The  diffraction  loss  term,  DFL,  is  the  same  term  used 

in  FACT  [5].   The  following  is  a  quote  from  Reference  5: 

A  preliminary  analysis  of  this  model  (Labianca,  1972) 
has  shown  that  the  leakage  term  does  not  correspond  to 
the  mechanism  which  dominates  duct  leakage  and,  in 
particular,  has  the  wrong  dependence  on  the  below-layer 
(thermocline)  gradient.   An  intensive  investigation  of 
all  available  surface  duct  models  has  been  completed 
recently  by.NUC  (Naval  Undersea  Center)  for  LRAPP  (Long 
Range  Acoustic  Propagation  Project) (Morris ,  1974).   This 
study  compared  several  models  with  each  other  and  the 
SUDS  (Surface  Duct  Sonar)  transmission  loss  data 
(Cummins,  1972).   Since  none  of  the  candidate  models 
emerged  as  clearly  superior,  AESD  (Acoustic  Environmental 
Support  Detachment)  is  currently  studying  possible 
modifications  to  the  surface-duct  equations.  [5] 
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Measured  data  taken  from  the  North  Atlantic  (ThorD  end 
Bernier,  1959  and  1962;  Sussman,  MacQonald,  and    Kanabis, 

1963),  Mediterranean  ( I.  al  lenient  and  Waterman,  1963;  Leroy, 
1964),  and  Northern  Pacific  (Greer  and  Rolam,  1954; 
Sheehy  and  Hal  ley,  1957)  f 3} . 


Thorp's  Ea 


Attenuation  = 


0.  IF1 


4  0F'~ 


1+40F^    4100+F'J 


x  10"°  dB/yd 


Fi  gure  16 
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11.  Convergence  zone  parameters  are  determined  at  the 
surface.   Errors  exist  when  source  and/or  receiver  are  at 
appreciable  depths. 

12.  When  CZ  conditions  are  not  present,  the  model 
assumes  there  is  no  deep  sound  channel. 

13.  Leakage  of  sound  into  a  shadow  zone  is  not  con- 
sidered.  If  there  are  no  CZ  conditions,  then  the  only 
illumination  of  the  shadow  zone  comes  from  "bottom  bounce. 
If  CZ  conditions  exist  the  shadow  zone  is  illuminated  by 

a  path  having  spherical  spreading  in  addition  to  the  bottom 
bounce  illumination. 
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V.   PROBABILITY  OF  DETECTION  (PD)  GENERATOR 

A.  INPUTS 

The  sub-program  for  probability  of  detection  contains 
values  from  the  cumulative  distribution  function  (CDF)  for 
a  standardized  normal  random  variable,  CDF  =  F(x).   Values 
of  x  corresponding  to  F(x)  =  0.5,  0.55,  ...,  0.95  are  stored. 

The  figure  of  merit  (FOM)  and  system  sigma  are  entered 
by  the  operator  as  part  of  the  transmission  loss  initializa- 
tion procedure.   System  sigma  is  defined  as  the  square  root 
of  the  sum  of  the  variances  of  the  individual  terms  in  the 
passive  sonar  equation: 

System  sigma  =  (a^  +  a^  +  a^  +  a^  +  a^ 

where  a  is  the  deviation. 

Rule-of-thumb  values  for  system  sigma  are  included  in 
Operator  Note  3  (Appendix  C)  for  those  cases  in  which  variances 
of  the  individual  term  are  unknown. 

B.  MODEL 

The  detection  variable  X  is  defined  as  the  ratio  of  the 
signal  excess, 

SE  =    FOM  -  TL  , 

to  the  system  sigma: 

X  =  SE 

System  Sigma 
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The  computed  X  is  compared  to  the  stored  values  of  x 
to  determine  probability  of  detection  at  which  PD  =  F(X) . 

C.   OUTPUT 

The  program  plots  a  probability  of  detection  (PD)  versus 
range  directly  below  the  TL  versus  range  plot  (Fig.  15). 
The  PD  routine  computes  PD  in  increments  of  0.05  so  that  a 
smooth  curve  does  not  result. 


VI.   CONCLUDING  REMARKS 

Advanced  development  of  an  on-board  environmental 
prediction  system  resulted  in  the  Integrated  Command  ASW 
Prediction  System  (ICAPS)  computer  program.   ICAPS  uses  the 
FACT  transmission  loss  program  so  that  it  requires  computer 
resources  comparable  to  those  available  on  aircraft  carriers. 
This  development  effort  is  continuing  with  emphasis  upon  the 
development  of  a  smaller  version  of  FACT.   A  viable  alternative 
would  be  an  augmented  version  of  PEAPS. 

Fleet  introduction  of  an  environmental  prediction  system 
need  not  be  delayed  until  a  foolproof  system  is  available. 
With  PEAPS  available,  even  in  its  present  limited  form, 
operational  evaluation  of  an  interim  system  with  some  known 
shortcomings  is  entirely  feasible.   A  continuing  dialogue, 
carried  on  between  the  developmental  agency  and  the  operating 
forces,  would  facilitate  implementation  of  system  improvements 
based  upon  fleet  experience. 

The  operating  forces  have  never  had  a  tactical  environ- 
mental prediction  capability  available  at  the  individual 
platform  level.   PEAPS  provides  an  inexpensive  way  to  evaluate 
the  operational  utility  of  such  a  system. 
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APPENDIX  A 
PEAPS  HARDWARE  DESCRIPTION 

For  this  program  the  Hewlett-Packard  Model  9  810A 
calculator  must  be  configured  with  the  following  options: 
Option  001:   Total  of  111  data-storage  registers; 
Option  003:   Total  of  2036  program  steps; 
Option  004:   Printer. 

Read-only-memory  units  (ROM's)  required  are: 
Mathematics  ROM  -  hp  -  112 10A 
Printer  Alpha  ROM  -  hp  -  11211A 
Plotter  ROM  -  hp  -  112 15A. 

The  Hewlett-Packard  Model  9  862A  Plotter  is  required  to 
produce  the  graphical  output. 
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APPENDIX  B 

The  following  pages  are  samples  of  the  Operator  Worksheet 
developed  for  use  with  the  PEAPS  program. 
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PEAPS  OPERATOR  WORKSHEET 

(See  PEAPS  Operator  Note  Sheets) 
BT  Information  DTG: 


Location: N/S , E/W 

\ 

1.  Turn  the  calculator  and  plotter  on.   "LOAD"  the  SSP/Ray 
trace  program  cards.   Put  the  plotting  paper  on  the 
plotter  and  set  the  "GRAPH  LIMITS"  of  the  plotter. 

2.  Press  "STOP,  END,  CONTINUE"  (this  sequence  is  used  when 
a  new  SSP  is  to  be  generated. 

Note  1:   If  you  enter  an  erroneous  value: 

a.  If  you  have  not  pressed  "CONTINUE",  press 
"CLR  X",  then  enter  the  correct  value. 

b.  If  you  have  pressed  "CONTINUE",  you  must 
re-enter  all  the  points  again;  press  "GO  TO, 
LABEL,  A",  then  continue  with  Step  3  below. 

3.  Press  "CONTINUE". 

4.  SS?  inputs:   ENTER  number  of  ET  points  (15  maximum) 

Press  "CONTINUE"  after  each  value  is  entered  to  advance 
the  program  _       0ut?ut 

input       Input       Input  Speed  of  Sound 

Point  Depth  (ft)   Temp  (°F)   Salinity  (ppt)   (ft/sec) 

1     Ch_0 

2 


10 


15 
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Note  2:   Change  the  plotting  paper  on  the  plotter.   Enter 
the  ray  trace  program  MODE  indicator. 

0   =  Automatic  compute  and  plot 

MODE:  1  -  Automatic  compute,  no  plot 

2  =  Manual  inputs  for  plotting,  Automatic 
compute . 

If  a  ray  trace  input  error  is  made,  press  "GOTO, 

LABEL,  B,  CONTINUE",  then  input  the  MODEL  indicator 

and  ray  trace  data  again.   (This  sequence  is  used 

to  generate  a  new  trace  using  the  old  SSP  data) 

Press  "CONTINUE" 

4.  Ray  trace  inputs: 

Press  "CONTINUE"  after  each  input  to  advance  the  program, 

Source  Depth: (ft) 

If  Mode  2  is  used: 

Minimum  Ray  Angle  x  2:  (degrees) 

Angle  Increment  x  2:    (degrees) 

Maximum  Ray  Angle  x  2:  (degrees) 

Note  3 :   If,  while  in  Mode  0   and  after  the  bottom  limiting 
ray  has  been  drawn,  a  channelled  ray  is  drawn  to 
the  plot's  range  limit;  see  Operator  Note  sheet 
section  Note  2E . 

Change  the  plotting  paper  on  the  plotter.   Press 
"STOP,  END".   LOAD  the  TLOSS/PD  program  cards. 
Press  "STOP, END".   (The  "STOP,  END"  sequence  is 
used  when  a  TLOSS/PD  input  error  is  made,  or  a  new 
TLOSS/PD  plot  is  to  be  generated  based  on  the  old 
SSP/Ray  trace  data.) 

5.  TLOSS/PD  inputs: 

Press  "CONTINUE"  after  each  input  to  advance  the  program, 

Press  "CONTINUE" 

Source  Frequency:  (Hz) 

Figure  of  Merit  (FOM)  : (dB) 

System  Sigma: (use  10  if  unknown) 

Wave  Height:  0   1  2  3  4  5  6  7  0  9  10  11  12  (ft) 

Receiver  Depth:        (ft) 

Bottom  Type:    123456789 


64 


Note  4 :   If  a  new  ray  trace  is  to  be  generated  using  the 
old  SVP  data  after  the  TLOSS/PD  program  has  been 
run,  Press  "STOP,  END".   LOAD  the  SVP/Ray  trace 
program  cards. 

Press  "GOTO,  LABEL,  B,  CONTINUE",  then  see  Note  2 
above . 
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APPENDIX  C 
PEAPS  OPERATOR  NOTE  SHEETS  FOR  THE  HP9  810A 

SSP/RAY  TRACE  PROGRAM 

The  SSP/RAY  Trace  program  is  recorded  on  two  short 
magnetic  cards,  each  card  having  two  sides.   The  program 
is  inserted  into  the  calculator  as  follows  (it  is  assumed 
the  calculator  is  turned  on  and  paper  has  been  placed  on 
the  plotter) . 

1.  Press:   END,  LOAD. 

2.  Insert  the  cards  starting  with  side  1  sequentially 
through  side  4. 

3.  Press:   END,  CONTINUE. 

After  the  SSP  axes  have  been  drawn  the  printer  will  indicate 
"SEE  NOTE  1",  referring  the  operator  ot  Note  1. 

Note  1 

A.   The  sound  speed  profile  (SSP)  generator  will  accept 
up  to  fifteen  (15)  points.   These  points  are  entered  from 
the  PEAPS  Computer  Operation  V7orksheet  starting  with  the 
first  point  being  the  surface  (depth  =  0)  temperature  and 
salinity.   The  data  are  entered  as  follows: 
1.   Press:   CONTINUE 

The  printer  will  request  "NO.  PTS".   Enter  the 
number  of  BT  points  to  be  entered  as  indicated  on  the 
work sheet. 
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2.  Press:   CONTINUE. 

The  printer  will  request  DEPTH.  Enter  the  depth 
in  feet  of  the  point  for  which  data  are  being  entered  (enter 
0  for  point  1) . 

3.  Press:   CONTINUE. 

The  printer  will  print  the  depth  entered  and 
request  TEMP.  Enter  the  temperature  in  °F  for  the  point 
of  interest. 

4.  Press:   CONTINUE. 

The  printer  will  print  the  temperature  entered 
and  request  SALIN.   Enter  the  salinity  in  parts  per  thousand 
for  the  point  of  interest. 

5.  Press:   CONTINUE. 

The  printer  will  indicate  the  salinity  entered 
and  the  speed  of  sound  in  feet  per  second  corresponding  to 
the  depth,  temperature,  and  salinity  values  immediately 
preceding  the  sound  speed  value.   The  plotter  will  plot  the 
speed  of  sound  indicated  on  the  printer  at  the  corresponding 
depth.   If  the  data  entered  are  from  a  point  less  than  the 
total  number  of  points  to  be  entered,  the  printer  will  indi- 
cate DEPTH.   Enter  the  depth  of  the  next  point  of  interest 
and  continue  from  step  3  above.   If  the  date  entered  is  for 
the  last  point,  the  program  will  discontinue  plotting  and 
the  printer  will  indicate  "SEE  NOTE  2".   This  refers  to 
Operator  Note  2  below. 

B.   If  the  BT  trace  does  not  reach  the  bottom: 
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1.  Use  at  least  two  points  from  historical  data 
and  enter  these  as  the  last  two  points.   The  more  points 
entered,  the  more  accurate  the  plot  if  the  bottom  is  con- 
siderably deeper  than  the  BT  trace  depth.   See  Section  G 
below  for  merging  data. 

2.  If  historical  data  are  not  available,  enter  the 
last  BT  temperature  as  the  bottom  temperature,  the  bottom 
depth  as  measured  or  approximated,  and  35  ppt  as  salinity. 

C.  If  salinity  values  are  not  available  from  measurement! 
or  historical  data,  enter  35  ppt. 

D.  If  you  enter  an  erroneous  value: 

1.  If  you  have  not  pressed  CONTINUE,  press  "CLR  X", 
then  enter  the  correct  value . 

2.  If  you  have  pressed  CONTINUE,  you  must  reenter 
all  the  points  again;  press  "GO  TO,  LABEL,  A",  then  continue 
from  step  2  in  Section  A  above. 

E.  If  you  have  finished  with  one  sound  speed  profile 
and  want  to  enter  a  new  one,  press  "END",  then  go  to  step  1 
in  Section  A  above. 

F.  Pressing  "CONTINUE"  after  each  input  advances  the 
program  and  enters  the  input  point  into  storage. 

G.  If  the  in-situ  and  historical  data  do  not  match 
closely  together,  merging  becomes  necessary.   There  is  no 
single  correct  method  for  merging  historical  and  in-situ 
data.   Suggested  procedures  are  to  fare  (draw  a  smooth  curve) 
the  two  profiles  together,  to  extrapolate  the  in-situ  trace 
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down  to  the  historical  data,  or  to  extend  the  profiles  in 
the  area  of  interest  until  they  meet.   Fig.  17a  shows  a  plot 
of  in-situ  SSP  and  historical  data  on  the  same  axis.   It  is 
suggested  that  this  be  done  for  all  SSP  plots  if  historical 
data  are  available.   A  trace  of  the  historical  data  need 
not  start  at  the  surface.   A  depth  500  feet  above  the  last 
BT  depth  should  suffice.   Figures  17b  and  17c  show  examples 
of  extending  the  two  traces  to  achieve  a  merge.   In  any 
event,  a  smooth  blending  of  the  traces  should  be  attempted. 
In  some  instances  it  may  be  necessary  to  try  more  than  one 
technique  in  order  to  achieve  a  reasonably  correct  facsimile 
of  the  sound  speed  structure. 

Note  2 

This  Note  deals  with  procedures  for  the  ray  trace 
generator. 

A.  Change  the  plotting  paper.  ; 

B.  Enter  the  ray  trace  program  Mode  indicator  (0,  1, 
or  2)  . 

Mode  0  gives  an  automatic  search  for  and  plotting 
of  convergence  zone  critical  rays.   The  minimum  source 
emission  angle  that  will  produce  a  bottom  bounce,  and  the 
minimum  (inside)  and  maximum  (outside)  convergence  zone 
annulus  ranges  (if  CZ  is  present)  are  stored. 

Mode  1  gives  the  same  computation  as  Mode  0  except 
that  the  ray  traces  are  not  plotted. 
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Mode  2 :   the  operator  enters  twice  the  minimum  and 
maximum  emission  angles  to  be  traced  from  the  source  along 
with  twice  the  incremental  angle  between.   Doubling  the 
desired  angles  and  angle  increment  is  necessary  since  the 
program  internally  halves  all  input  angles  prior  to  conver- 
sion to  radians  and  storing.   When  the  manual  input  plotting 
is  completed,  the  program  changes  to  Mode  1  to  determine 
the  critical  parameters  to  be  used  by  follow-on  programs. 

C.  Program  Mode  0  is  intended  to  be  the  primary  Mode 
for  fleet  use.   Mode  1  is  intended  for  use  when  a  ray  trace 
is  not  desired.   Mode  2  is  intended  to  be  used  for  ray  trace 
plots  of  specific  source  emission  angles  as  well  as  for 
determination  of  critical  parameters  for  follow-on  programs  ■ 
it  is  foreseen  as  a  training  aid. 

D.  If  you  have  made  an  incorrect  input  and  entered  it 
with  the  CONTINUE,  press:   GO  TO,  LABEL,  B,  CONTINUE. 
Re-enter  the  ray  trace  input  parameters  starting  with  the 
Mode  indicator. 

E.  It  is  possible,  under  certain  sound  velocity  profile 
conditions,  for  the  ray  trace  program  in  the  automatic  modes 
to  trace  a  ducted  ray  to  the  end  of  the  plot  while  searching 
for  the  convergence  zone  critical  ranges.  If  this  situation 
should  occur  you  must  check  the  stored  convergence  zone 
ranges  as  follows: 

1.   Press:   X«- (  )  ,  008,  t,  2000,  r. 
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This  calls  out  the  minimum  convergence  zone 
range  and  converts  it  from  yards  to  nautical  miles. 

2.  Compare  the  value  in  the  Y-register  with  the 
minimum  convergence  zone  range  indicated  on  the  ray  trace 
plot.   If  they  differ  considerably, 

a.  Enter  the  minimum  convergence  zone  range 
from  the  ray  trace  plot; 

b.  Press:   t,  2000,  x  (multiply),  Y->  (  ),  008. 
This  converts  the  ray  trace  plot  minimum 

CZ  range  from  nautical  miles  to  yards  and  stores  it  in  the 
appropriate  storage  register. 

3.  Press  X«-  (  ),  009,  t,  2000,  t. 

This  calls  out  the  maximum  CZ  range  and  converts 
it  to  nautical  miles. 

4.  Compare  the  value  in  the  Y-register  with  the 
maximum  CZ  range  off  the  ray  trace  plot.  If  they  differ 
considerably; 

a.  Enter  the  maximum  CZ  range  from  the  ray 
trace  plot. 

b,  Press:   i,  2000,  X  (multiply),  Y+ (  ),  009. 
This  converts  the  ray  trace  plot  maximum 

CZ  range  from  nautical  miles  to  yards  and  stores  it  in  the 
appropriate  storage  register. 

5.  Proceed  with  the  transmission  loss/probability 
of  detection  program. 
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F.  If  you  want  to  do  another  ray  trace  and  plot  using 
the  same  sound  speed  profile  but  different  input  parameters: 

1.  Press:   GO  TO,  LABEL,  B,  CONTINUE. 

2.  Enter  the  new  ray  trace  inputs  starting  with  the 
Mode  indicator. 

G.  "SEE  NOTE  3"  and  "END"  will  be  printed  on  the  output 
tape  and  the  calculator's  X,  Y,  and  Z  registers  will  illuminate 
when  the  ray  trace  portion  of  the  program  is  completed. 

H.   "NO  CZ",  "SEE  NOTE  3",  and  "END"  will  be  printed  on 
the  output  tape  to  indicate  that  the  program  has  ended  and 
there  is  no  convergence  zone  present. 

I.   If  "NO  CZ",  "SEE  NOTE  3",  and  "END"  is  printed  out 
and  you  would  like  a  ray  trace  plot,  you  will  have  to  use 
Mode  2  and  manual  emission  angle  parameters  to  get  a  plot 
(see  E  above).   At  the  end  of  the  Mode  2  plots  "NO  CZ", 
"SEE  NOTE  3",  and  "END"  will  be  printed  on  the  output  tape. 

J.   If  another  sound  velocity  profile  is  to  be  entered 
for  subsequent  ray  trace  plots,  press:   END,  CONTINUE,  to 
start  the  entire  sequence  again. 
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TRANSMISSION    LOSS/PROBABILITY    OF    DETECTION    PROGRAM 

Note    3 

A.  Change  the  plotting  paper. 

B.  The  TL/PD  magnetic  cards  are  loaded  in  the  same 
manner  as  the  SSP/Ray  Trace  cards. 

C.  Press:   STOP,  END  to  initialize  the  program  function. 
This  sequence  is  used  when  a  TLOSS/PD  entering  error  is  made, 
or  a  new  TLOSS/PD  plot  is  to  be  generated  based  on  the  same 
SSP/Ray  Trace  generated  data. 

D.  Source  frequency  is  to  be  entered  in  Hertz  (Hz) . 

E.  Figure  of  Merit:   ensure  that  you  use  all  the  same 
units  of  reference  (yPa,  or  yb) .   DO  NOT  add  100  to  the 
computed  FOM. 

F.  System  Sigma 

If  the  variances  of  the  sonar  equation  terms  are 
known,  enter  the  square  root  of  the  sum  of  the  variances. 
If  the  variances  are  unknown,  choose  one  of  the  following: 

1.  Use  6  dB  if  ambient  noise  measurements  have  been 
made,  and  submarine  speed  and  type  are  known. 

2.  Use  8  dB  if  ambient  noise  is  estimated  from  a 
forecast,  submarine  speed  is  known  to  within  3  kts ,  and 
submarine  type  is  known. 

3.  Use  10  dB  if  ambient  noise  is  estimated  from 
forecasts  and  submarine  speed  and  type  are  uncertain. 
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G.   Wave  Height  (IWH) 


Sea  State 

Term 

He. 

ight  of  Wave 

(ft) 

1 

Calm 

0 

to 

1 

2 

Wavelets 

1 

to 

2 

3 

Slight 
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to 

4 

4 

Moderate 

4 

to 

8 

5 

Rough 

8 

to 

13 

6 

Very  rough 

13 

to 

20 

7 

High 

20 

to 

30 

8 

Very  high 

30 

to 

45 

9 

Phenomenal 

over  4  5 

H.   A  bottom  type  indicator  (MGS)  of  1  indicates  little 
sound  loss;  a  bottom  type  indicator  of  3  indicates  a  large 
sound  loss. 

I.   If  new  TLOSS/PD  data  is  to  be  entered  to  make 
additional  plots  from  the  same  ray  trace  data,  press  STOP, 
END.   Change  the  plotting  paper,  then  press  CONTINUE.   Then 
enter  the  new  data. 

J.   If  a  new  ray  trace  is  to  be  generated  from  the  old 
SSP  information  after  the  transmission  loss/probability  of 
detection  program  has  been  run:   press:   STOP,  END,  LOAD. 
Load  the  SSP/Ray  Trace  program  using  the  magnetic  cards, 
change  the  plotting  paper,  press:   GO  TO,  LABEL,  B,  CONTINUE, 
then  follow  the  instructions  in  Note  2. 
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APPENDIX  D 

Flow  charts,  program  listings,  and  output  from  the 
HP9810A  and  IBM  360  are  included  within  this  Appendix. 
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Speed  of  Sound  Profile  /  Ray  Trace  Program 
flow  charts  and  sample  outputs 
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(/*)  SSP  critical  points 
selection  program 
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NG  =  NG  +  1 
increment 
storage 
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(*0  SSP  critical  points  selection  program 


< 


1 

No 

store  Sy 

increment 

storage 

NG  =  NG  +  1 

set  the  FLAG 

1 

' 

& 


store  Sx 
set  the  FLAG 


<} 


store  Sy 

increment 
storage 

NG  =  NG  +  1 


& 


Run  = 
Run  +  1 


Sy  -  Sx 
Sx  =  the 
next  S 


(*0  SSP  critical  points  selection  program 


0 


clear 
the  FLAG 


store  Sy 

increment 
storage 

NG  -  NG  +  i 


Sy  =  Sx 
Sx  =  the 
next  S 


Yes 


store  Sy 

increment 
storage 

NG  =  NG  +  1 


<} 


store  Dl  and  SI 
in  GDI  and  CS1 

store  the  last 
D  and  S  in  CD 6 
and  GS6 


I 


LABEL  B 
see  note  2 


O 
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A  )   Speed  of  Sound  Profile  /  Ray  Trace 
Program 


s^>0 
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P.. 

J 


B. 

Fl 

n 

Z 

Z 
u 

U 
EI 
Z 

E      - 
C    B. 


w 

a 

o 
o 

{-I 

E-» 

& 


*A 


B.. 


t  I  !  i  i- 


+-H_4-M_4HH-4_4-4~^4H^^ 


b1 


14 
IS 


B 

1:1 


taai  x  .Ln  hms<3 
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(6)      Convergence  Zone  Test 


m  =  2 


No  CZ 


Compute 
SV,    THP, 
THPP 


Compute 
SV,    THP, 
!£HPP 


<} 


m  ■  m+1 


<} 


Yes 


<H? 


es 


87 


(7)   Computation  of  SV ,    THP,    THPP 


SV   =    (SD-Dm-l)x 
|.Sra-Sm-l\     , 
I  Dm-Dm-l] 

Sm-1 


IHP  - 

2SV 
CS2 

f 

THPP  = 

V2~ 

2SV 
CS6 

^ 
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Speed  of  Sound  Profile  /  Ray  Trace 
Program 


Ahmin=thpp\ 

/set  plotting\ 


o 


Trace  t he- 
Ray 


THMIN=THMIN+.5° 
set  plotting 
indicator 


9?HMIN=$HMIN-.l 

clear  plotting 
indicator 


set  FLAG 


?EBB=!rKMirr 
set   ICZ 
set  plotting 
indicator 
CZmax-R 


ie  ray 


Draw  t 
that  just 
grazes  the 
bottom 


THMIN  - 
THP,  set 
plotting 
indica- 
tor 
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(8)  Ray  Trace  Program 


No 


TH  = 

THMII' 

z  -  sd 

R  =  ff 


Clear 
Plotting 
Indicator 


A  Surface  Re- 
flected Ray, 
Return  to  Call 
Point 


RADINV 


n2 


-G(l  -  (TKc/2j) 

fG(Z  -  D(m-l))| 
+  S(m-l) 


DR  = 


5^>(m)-D(m-l) 


A  Bottom  Reflected 
Ray,  Return  to 
Call  Point 


90 


(8)  Ray  Trace  Program 


© 


RDR  = 
RADINV x 
DR 


I 


R  = 
R+DR 


Z  = 


Z+S|R  +THDR 


Trace  the 

Ray 
Segment 


A  Channeled 
Ray,  Return 
to  Call  Point 
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J 
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14 


KMI   X  U3  BlcBfl 
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T 


n. 
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K 
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B 
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cd 
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& 

rH 
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(0 

cd 
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Mill- 
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id 

N 


Id 


B 
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11 


CEHI  X  Uj  HU13<I 
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0 


D  )   Speed  of  Sound  Profile  /  Ray  Trace 
Program 


O 


1 


Trace  the 
Ray 


CZMIN  =  R 

THMIN=THMIN+.5° 

set  plotting 
indicator 


i. 


Trace  the 
Ray 


(U) 


es 


£> 


l'HMIN=TKMIN+.l 

clear  plotting 
indicator 


CZMAX 
R 
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P. 

e  >> 
w  u 


1-4^.. ..{....I  ;   )   ;  ^HHHg^H^^Hm^.;^^.H^f-4HH"-HH 


LI 


M 


LE05I   X  JU3  RLrDC! 


96 
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0  Speed  of  Sound  Profile  /  Ray  Trace 
Program 


THMIN=THMIN-.4( 

clear  plotting 
indicator 


<} 


Trace  the 
Ray 


es 


CZMIII  =  R 
rHMIN=THMIN+.l° 


•o 


THMIN=THMIN-. 

set  plotting 
indicator 


Trace  the 
minimum  CZ 
range  ray 


(12) 


9R 


© 


C  )  Speed  of  Sound.  Profile/  f  F 
Ray  Trace  Program 


Trace  the 
Ray 
(MODE  2) 


THMIN  = 
THMIN  +  DTH 


No 


is 

iIN>TJ 
? 
'Yes 


set 

M0DE=1 


set  IBB 

convert  CZMIN 
and  C2MAX  froir 
feet  to  yards 


set  the 

plotting 

indicator 


See  Note  3 


(13) 
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PEAPS  HP  981 OA 
Speed  of  Sound  Profile/Ray  Trace 
Program  Listing 


101 


0  0  0  0  -  - 

GTO- 

--44 

0001 -- 

4.  - 

--04 

0002-- 

cr 

.j       — 

--05 

0  O  0  3  -  - 

YTO- 

--40 

0004-- 

IND- 

--31 

0  0  0  5  ~  ~ 

k     - 

-  -  1  4 

0  0  06~- 

3    - 

-  -03 

0  U  07-- 

xfu- 

—  _  'T'  O 

0008-- 

-    - 

--34 

0  O  O  9  ~  - 

o.    - 

-  -  1  3 

0010  — 

6    - 

-  -  0  6 

0011  — 

XTQ- 

, IM  .„  '"i  *'_i 

001 2-- 

-    - 

--34 

0  0 13  — 

b    - 

--14 

0  014  — 

XFR- 

--67 

0  0  1  5  -  - 

IHH- 

--31 

O  0 16  — 

o.    — 

--13 

0017  — 

X  T  d  - 

*~~  """  il~  O 

0018  — 

INH" 

-  -  3  i 

0019  — 

b     - 

—  14 

0  0  2  0  -  - 

3    - 

-  -  O  3 

0021  — 

XTO- 

—  —  ■*•  o 

0022--- 

+    - 

—  3  3 

0023-- 

ij.    — 

--13 

0  0  2  4  -  - 

i     — 

-  -  0  7 

0025-- 

XTO- 

ti,.  o 

0026--- 

+    - 

_  ~    -t  ^ -( 

0  0  2  7  -  - 

b     - 

--14 

0  0  2  8  -  - 

1     - 

-  ~  0  1 

0  0  29-- 

XTQ- 

■*"*"£,  O 

0  0  3  0  — 

+    - 

-  -  3  3 

0  0  3  1  -  - 

2    - 

--02 

O  0  3  2  -  - 

S/R- 

{  : 

0  0  3  3  -  - 

XFR- 

—  67 

0  0  3  4  -  - 

IND- 

-  -  3  1 

0  035--- 

o.    — 

-  - 1  3 

0  0  3  6  -  - 

UP- 

—  2  i* 

0  0  3  7  — 

2     — 

--02 

0  0  3  8  -  - 

X  7  0  - 

L—    •-.' 

'0039-- 

+ 

-  -  3  3 

0  0  4  0  ~  - 

ft      — 

—  1 ::: 

0  0  4  1  -  - 

XFR- 

—  6  7 

0042-- 

IND- 

-  -  3  1 

0  0  4  3-- 

o.    — 

--13 

0  0  4  4  -  - 

S--R- 

r  ,' 

O045-- 

FMT- 

-  ~  4  2 

0  0  46  — 

1    - 

--01 

0047-- 

8    - 

-  -  j  y 

0048-- 

0    - 

-  ~  0  0 

■  y  y  4  9  -  - 

cr 

.j      — 

-  -  y  5 

0050.-- 

2 

--02 

0  0  5  1  -  - 

1     - 

-  -  0  1 

0  U  ■:•  £  — 

0      ~ 

-  -  0  0 

0053  — 

UP- 

""if 

0O54---  4 0  4 

■  6 06 

0 00 

0   0  0 

FMT 42 

1 01 

2  02 

5 05 

-  -  0  0 

--00 


0O55- 
0  0  56- 
0  057- 
0  0  5  8  - 
0  0  5  9  - 
0  0  6  0  - 
0  0  6  1  - 
0062- 
0  0  6  3  - 
0  0  64- 
0  0  6  5  - 
0066  - 
0067- 
0  0  68- 
O069- 
0  0  70- 
0071- 
0072- 
00  73- 
0  0  ?  4  - 
0  0  75- 
00  76- 
0  0  7  7  - 
0  0  78- 
0  079- 
0080- 
O081  - 
0082- 
0  0  8  3  ~ 
0  0  84- 
0  0  8  5  - 
0  0  8  6  - 
00  37- 
0  0  8  8  - 
0  0  8  9  - 

0  090- 

0  0  9  1  - 
00  92- 
0093- 
0  0  9  4  - 
0  0  9  5  - 
0  0  9  6  - 
0  0  9  7  - 
009  8  - 
0  0  9  9  - 
0  1  0  0-  - 

0101- 

0102- 

0  1  0  3  ~ 
0104  - 
0  10  5- 
0106- 

0107  - 


y 

0 

UP 

HS' 


—  02 


-EEX- 
_  4  

-FMT 

-  0   

-  UP 


-  0   - 

-  0   - 

-  F  M  T  - 
_  i    _ 

-  UP- 

-  1  - 

-  0   - 

-FMT- 

-  1  - 


-  0   

-  0   

-FMT--. 

-  6  — 

-  0 

-FMT--' 

-  1  --■ 

-  1  --• 

-  1  --■ 


UP 

4 

9 

0 
M 


04 
4  2 

01 

03 

00 

'  iZ.    i 

•04 

-  0  0 
"00 

-42 
-01 

■01 
■00 

•42 

■01 

•05 

0  5 

■00 
■00 

■42 

■01 

06 

•21 

01 
00 

■42 

■01 

01 
•01 

•02 
0  7 
05 

04 

1 1 

0U 
0  0 


0 

1  0  8  - 

-FMT- 

--4  2 

0 

[09- 

-     1     - 

--81 

0 

1  lO- 

-    UP- 

ti  i 

0 

ll  1- 

-FMT- 

—  t  c 

0 

1  12- 

-     1     - 

--01 

0 

1  13- 

-FMT- 

--42 

0 

l  14- 

-VT0- 

7-40 

0 

L  15- 

-    if    - 

--56 

0 

1.16- 

-    E    - 

--60 

0 

.  j.  i 

-    E    - 

-  -  6  0 

0 

[18- 

-  n  - 

-  -  6  3 

0 

[19- 

-CHT- 

-  -  4 ":;' 

0 

[20- 

- !::!  <  Y  - 

_.  _  cj  •'• 

0 

.   i™   J. 

-    F    - 

—  16 

0: 

2  2  — 

-XT0- 

—  —  2  '■' 

0 

2  o  — 

-IH  V  - 

-  -  3  5 

0 : 

[24- 

-YTU- 

--40 

6 : 

125- 

-    E    - 

-  -  6  0 

0 : 

.  2  6  - 

-  c  - 

~  -  6  i 

0: 

.  c.  f   — 

-CNT- 

0 ; 

2  8_ 

-    YE- 

--24 

0 : 

2  9  - 

-CNT- 

- 47 

0 : 

,  3  0  - 

-     1     - 

-  -  0  1 

O : 

[31- 

-      0      - 

-  -  0  0 

0  j 

•~i  '"i  .., 

-      0      - 

-  ■-  0  0 

0: 

.  •-'  •-'  "™ 

-PSE- 

-J  f 

0 : 

.34- 

-FNT- 

--42 

0 : 

.35- 

-    4    - 

-  -  0  4 

0 : 

.  3  6  - 

-    6    - 

-  -  0  6 

0 : 

.  •„"•  i''  ■• 

_      Q      _ 

~~l  i 

0: 

.  3  8  - 

-      0      - 

--00 

0  ' 

3  9  - 

-  :=•=:  t  ij  - 

£m    *—' 

0 : 

'■40- 

-    o.    - 

-  -  1  3 

0 : 

41- 

-    4    - 

-  -  0  4 

0] 

•i  -j  _ 

—    t     — 

-  ■■-  0  7 

0 : 

43- 

-  >',  7  d  ~ 

-  -  2  3 

0 : 

44- 

-    b     - 

--14 

0  j 

45- 

-    3    •- 

...  ..  y  3 

0 : 

4  6  - 

-      0      - 

-•--011 

0 : 

47- 

-      0      •- 

—  0  0 

0  J 

4  8  - 

-    UP- 

mm,    ....    '"I1    "7 

0 : 

49- 

—      ft 

-  - 1  3 

0  3 

50- 

-FMT- 

-  -  4  2 

01 

51- 

-    1    - 

--01 

0  j 

cr  •-■  _ 

-  IJ  p  - 

iL.  i 

0  J 

•j  3  — 

-    b    - 

--14 

0  J 

54- 

-FMT- 

--42 

0  J 

s  s  — 

-     1    - 

--01 

0  \ 

56- 

-PUT- 

...  _.  4  5 

0  \ 

57- 

—    1    — 

--01 

0  1 

cr  o 
•J  o  — 

-      0      - 

-  ~  0  0 

0  \ 

59- 

-      0      ~ 

■  ~  0  U 

0  \ 

6  0  - 

-kto- 

—  2  3 

0  1 

61- 

--    +    - 

_..  .—  ^  ^  -( 
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6162  —  o.  1  3 

0  2 1 6  -  -  6  0  6 

8278--LBL- 51 

0  1  6  3  -  -  1  01 

O  2  1  7  -  -  O  0  0 

8  2  7  1  -  -  H  62 

0  1  6  4  -  -  X  TO 2  3 

O  2  1  8  -  -  O  80 

0  2  7  2  -  -  C  L  R 28 

0  1  6  5  -  -  +  3  3 

0219  —  CH  3 3  2 

6273--  9  11 

0  1  6  6  -  -  b 14 

0  2  2  O  -  -  X  TO 2  3 

8274--  9 11 

01 67--  b  14 

O  2  2  1  -  -  o. 1  3 

8275--  UP 27 

01  68--  UP 27 

0222--  2 02 

8  2  7  6  -  -  C  L  X  -  -  -  3  7 

0  1  6  9  -  -  5  0  5 

0  2  2  3  -  -  5 O  5 

O  2  7  7  -  -  X  T  0 2  3 

0  1  7  0  -  -  3 0  3 

O  2  2  4  -  -  X  T  0 2  3 

0273--  I  ND 31 

0171  —  X>  Y 53 

0225--  b  14 

0  2  7  9  -  -  o.  1  3 

0  1  7  2  -  -  U 0  0 

0  2  2  6  -  -  o. 13 

0  2  8  0  -  -  1 61 

01 73--  1 01 

0227--  UP 27 

6231— XT  0 23 

0  1  7  4  -  -  4  04 

0228--  4 84 

8  2  8  2  -  -  +  33 

0175--  5 05 

O  2  2  9  -  -  6  0  6 

62  83  -  -  o.  1  3 

0176--FMT 42 

02:38--  5  05 

0284--  o.    1  3 

0177-  -  1  01 

0  2  31--  8  0  8 

8  2  8  5  -  -  X  <  Y 5  2 

01  78--  9  11 

0232  —  FMT  — -42 

8  2  8  6  -  -  8  68 

0  1  7  9  -  -  C  H  S 3  2 

0233--  1  0  1 

8237--  2  ---82 

0  i  8  0  -  -  i  —  0 1 

8234--  UP 27 

0283--  7  87 

0  1  8  1  -  -  2  82 

82  35-~  b  1 4 

8  2  8  9  -  -  6  8  6 

8182--EEX 26 

0236  — FMT 42 

0  2  9  O  -  -  C  L  P 2  8 

8  1  8  3  -  -  3  0  3 

823  7--  1 81 

8  2  9  1  -  -  F  M  T 4  2. 

01  84--  UP 27 

6238  — P  NT- 4  5 

82  92--  1  61 

0135  -  -  4  O  4 

823?--  2  - 8  2 

829  3  -  -  8  16 

0  1  8  6  -  -  6  0  6 

824  8--  5  85 

0294--  2 62 

0  i  8  7  -  -  3  0  3 

8  2  4  1  -  -  8  - 8  8 

0295  — FMT 42 

0  1  8  8  -  -  U  0  U 

82  4  2--  8  88 

6296--FMT 42 

0  1  8  9  -  -  F  M  T  -  -  -  4  2 

8  2  4  3  -  -  K  T  0 2  3 

0  2  9  7  -  -  Y '!'  0 4  6 

0  1  9  0  -  ~  1  8  1 

8244--  -  34 

6  2  9  8  -  -  Y  T  0 4  6 

0  1  9  1  -  -  U  P 2  7 

0245--  o.    1  3 

6  2  9  9  ~  -  fi  5  6 

0192--FMT- 42 

8  2  4  6  -  -  2  6  2 

0  3  O  O  -  -  C  L  F: 28 

0  1  9  3  -  -  1  0  t 

8  2  4  7  ~  -  5  8  5 

6  3  8  1  -  -  Y  TO 4  6 

0194--FNT 42 

0248  — XT  0 23 

6  3  6  2  -  -  E  66 

0  1  9  5  -  -  ii  6  3 

8  2  4  9  -  -  +  33 

8  3  W  3  -  -  E  fa  0 

0  1  9  6  -  -  E  6  U 

8258--  b  - -14 

8384  --CNT 4'7 

0  1  9  7  -  -  fi  56 

825 1--  b  14 

0  3  0  5  -  -  N  7  3 

01  98-- X  TO 2:]: 

8 2 5 2 - -  UP 2 7 

6  3  6  6  -  -  0 7  1 

0199  —  H  74 

8  2  5  3  -  -  2  62 

6  3  8  7  -  -  X  T  0 2  3 

O 2 O U  -  -  C  NT 4  7 

6254--EEX 26 

8308-^-  E  6  0 

0201--XCY 52 

8255--  2  82 

03O9--CNT 47 

0  2O2--  F  16 

6  2  56--X>Y- 53 

6  3  1  6  -  -  1  01 

0  2  0  3  -  -  X  T 1 3 2  3 
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Transmission  Loss  /  Probability  of  Detection  Program 
flow  charts  and  sample  outputs 
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TLOSS  /  PDET  Program 
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RS    = 
arccos 


CS(1 


Ml)  x 

CS(2)IA 
(2CS(l)  /  3)  X 
(CD(2)  -  CD( 1) 
(CS(2)   -  CS(1) 


RS    = 


CD(2) 


6  arccos  ^2) 


SSL    » 

.fifv/HxSFl     '?1 

ES"   iooo 


G   = 


CS(3)-CS(2) 

cLT3>rcDT27 


120 


(18)   Computation  of  RL2 
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(19)    Computation  of  RPT 
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I  )   TLOSS  /  PDET  Program 


Determine  the 
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factor  BT 
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(20)  Bottom  Loss  Factor  BT  Determination 
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(21 )  Determination  of 
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(22)  Computation  of  Absorption  Coefficient  (A) 
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(23)   Preset  ICZRC,  N,  DRWRK,  DRPERM 
Compute  THBB,  THBL 
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TLOSS  /  PDET  Program 
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(2*0      Bottom  Bounce   Path  TL03S  at  Range  R 
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(25a)   STL  m   STL  *  TL  Subprogram 
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(26)  Reliable  Acoustic  Path  TLOSS  at  Range  R 
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(2?)     Convergence  Zone  Path  TLOSS  at  Range  R 
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(28)  Deep  Sound  Channel  TLOSS  at  Range  R 
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TLOSS  /  PDET  Program 
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(30)  Determination  of  Probability  of  Detection 
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PEAPS  HP  981 OA 
Transmission  Loss/Probability  of  Detection 
Program  Listing 
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1  6  3  2  ~  - 

DN  — 

—  ■"■  s 

1  6  3  3  -  - 

+  -  - 

-  3  3 

1  6  3  4  -  - 

Y  T  Li  -  - 

-48 

1  6  3  5  -  - 

< 

-07 

i  6  3  6  -  - 

cr   _  ._ 

-05 

1  6  3  7  -  - 

G  T  0  -  - 

-  4  4 

1  6  3  8  -  - 

S/R  — 

r  i 

1  6  3  9  -  - 

1  -~ 

-01 

1  f,  4  g  -  - 

r   "™ 

-  0  7 

16  41  -  - 

c  L  x 

—  ■_'  i 

1  6  4  2  -  - 

up- 

"~  C  i 

1  6  4  3  -  - 

XFR  — 

...  f.f  (- 

1644  — 

•~i   tmm 

-02 

16  4  5  — 

6  — 

-0  6 

1  6  4  6  -  - 

v  =  y  -  - 

-50 

16,47-- 

1  -- 

-01 

1  6  4  8  -  - 

8  -  - 

-10 

1649--. 

v   _.  __ 

-02 

1650  — 

t 

-07 

1651  — 

XFR 

-  6  7 

1652-- 

{    ~  — 

-  0  7 

1653  — 

8  -  - 

-10 

1654  — 

UP  — 

— .  O  .•' 

1  6  5  5  -  - 

XFR— 

-  6  7 

1656  — 

{ 

-07 

1657-- 

i 

-07 

•j  r  cr  O 

x  

-  3  6 

1  6  5  9  -  - 

G  T  0  -  - 

-  4  4 

1  6  6  8  -  - 

S/R  — 

i  i 

1661  — 

3  -  - 

-  0  3 

1  6  6  2  — 

DN  — 

1  6  6  3  — 

+  -- 

...  o  ■— ' 

1664  — 

YTQ— 

-40 

1665  — 

"7   

-07 

166  6  — 

cr 

-05 

1  6  6  7  -  - 

XFR-- 

-  6  7 

1  6  6  8  — 

(' 

-07 

1669  — 

r   —  — 

-07 

1676  — 

UP  — 

_  *"'  "? 

£.  1 

1671  — 

XFR  — 

-67 

1  b  (   iC  —  — 

C|   

-1  1 

1  6  7  3  - 

3  -- 

-  0  3 

1674--; 

.•  **.  ■.  .* 

cr  ■-. 
—  ■_'  -I1 

1  £.  "3  C 

1  t'  (   -J 

1  -- 

-01 

1  6  7  6  -  - 

8  -- 

-10 

1  6  7  7  -  - 

2  -- 

-02 

1  6  7  8  -  - 

3  -- 

-03 

1679--: 

<FR  — 

-  6  7 

1  6  8  0  -  - 

y   

-1  1 

1681  — 

4  _„.. 

-0  4 

1  6  8  2  — ', 

,; ,:;  y  _  ... 

—  5  2 

1  6  8  3  — 

1  ~- 

-01 

1684  — 

r 

-07 

1  6  8  5  -  - 

9  — 

-1  1 

1  6  8  6  — 

8  -  - 

-10 

1  6  8  7  -  - : 

<:fr— 

-  6  7 

1  6  8  8  -  - 

-  --- 

-  3  4 

1  6  8  9  -  - 

9  -- 

-1  1 

1  6  9  0  ~  - 

3  -  - 

~  0  3 

1691--' 

,;  f  o  -  - 

—  2  3 

1  6  9  2  -  - 

ij.  -  - 

- 1  3 

1  6  9  3  -  - 

UP-- 

—  2  r 

1  6  9  4  -  - 

1   -- 

-01 

1  6  9  5  -  - 

0   -  - 

-0  0 

1  6  9  6  — 

[i  i  y  -  - 

-  2  5 

1  6  9  7  -  - 

r'T0--~ 

-  4  8 

1  6  9  8  -  - 

■Z|   

-1  1 

169  9  — 

6  -  - 

-0  6 

1  7  0  0  -  - 

:.ht— 

...  4  7 

1701  — 

1   -- 

-01 

1702  — 

IJ  F:'  -  - 

""  tl  f 

1703  — 

4FR  — 

-  6  7' 

1  7  0  4  -  - 

3  --  - 

-  0  3 

1  7  0  5  -  - 

■\  >  Y  -  - 

-  5  3 

1  7  0  6  ~  - 

1  -- 

-01 

1707  — 

\ 

-  0  7' 

1 7  0  8  -  - 

Cm 

-  0  2 

1  7  0  9  ~  - 

4  _- 

-04 

171  0  — 

:,;  p  R  _  ... 

-67 

171  1  — 

q  — 

-1  1 

1712  — 

•~i 

-03 

1713  — 

ao 

'"'  O 

1714  — 

i 

-07 

1  7  1  5  -  - 

"7   

-07 

1  7  1  6  -  - 

i'  T  0  -  - 

-40 

171 7-- 

+  -- 

i  r  i  y 

■-i 

-  0  3 

1719-- 

:  t  o  -  - 

-  4  4 

1720  — 

i  

~01 

1721  — 

0       m~.     _ 

-  0  3 

^  ~ 7  •"I  -"i 

J.  r  iZ.  ll 

4  -.. 

-0  4 

1  f  C.  -1' 

6  — 

-86 

1724-- 

:hs— 

mmt  **j  *"i 

1  ~7  •-!  cr 

1  1  l!  ■..' 

4  __ 

-04 

1  7  2  6  -  - 

,  -  - 

-  2  1 

lie!  r 

cr: 

-05 

i  "7  ■:>  o  „  . 
i  i  i^.  •_• 

1729-' 

1730-' 

1731- 

1  7  3  2  -  ■ 
1  {  3  3  — c 
1  7'  3  4  - 
1735- 
1  7  3  6  - 
1737- 
17  38- 
1733- 
1740- 
174  1- 


EEX' 

CHS' 

4  ' 
UP 


dn- 

l'/X' 


--84- 

Cm    t 

--01 

— 0  5 


1742--  UP 2 


A   Z> 


1744 

1745 

1746' 

1747 

1  7  4  8 ' 

1  7'  4  9 

17  5  0 

1751 

1752 

1753 


—  i~, 


Y  T  0 

b  ■ 

'...'  r-  ri 
< 

UP 

XFR 


UN' 
XEY' 

r   ■ 


1  7  5  4  -  -  9 

1755 3 

1  7  5  6  -  - 

1757- 

1758- 

1759' 

1  7  6  0  • 

1761- 

1762' 

1763' 

1764' 

1765' 

1  7  6  6  ■ 

1  f  b  r  • 

1768' 

1769' 

1770' 

1771' 

1  i  (  c ' 

17  73 

1774' 

1775 

1776' 

17  77 

1773 

1779 

1780' 

1  7  8  1 


2i 
0  3 
0  6 
3  6 
48 
14 
6  7 
■07 
•0  7 


1  1 


-i  t- 
-  3  0 


.•■  b 


1 


DIV 

1 

■DIV 
DN 

XEY 
UP 
DN 

•DIV 
■  DN 

RUP 
RUP 
XEY 
CHS' 

DIV' 

•  DN- 

•  .1 


01 

10 


-i  c.- 
C   J 


30 


<-i  cr 


:-:m 


0  2 
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•\    -?  I*i  •-! 
1  I  •->•:- 

X 36                   1 

1  "7  O  O '• 

'  T  0 4  0                   1 

1  7  8  4  -  - 

-    34                   1 

1785  — 

7    0  7                   1 

1  7  8  6  -  - 

5    0  5                   1 

1  707 t 

;  N  T 4  7                 1 

1  "7  O  C'  _  _ 

1   ---01            1 

1789--) 

•i  i  0  —  —  2  3         x 

1790 - - 

+ 3  3                   1 

1791  — 

9    11                   1 

1792  — 

8 10                   1 

1  7  9  3  -  - 1 

-T0 44                   i 

1 7  9  4  — 

1 01            1 

1795  — 

3    10                   1 

1  7  9  6  -  - 

2    02                  1 

1  7  9  7  — 

2    02                  1 

1  798  — 

1    -01                   1 

1  7  9  9  — ! 

<T0 23                   1 

1  8  0  8  -  - 

9 i  i                   1 

1  '3  01  — 

8    10                   1 

180  2--! 

no 23                 1 

1  8  0  3  -  - 

3    0  3                  1 

1  8  0  4  -  - 1 

JNT---47                  1 

1  3  0  5  -  - : 

<FR -67                  1 

1  8  0  6  -  - 

9    11                   1 

1  8  0  7  -  - 

7    07                  1 

1  8  0  8  -  - ! 

•no -23                 1 

1  8  0  9  -  - 

9    11                   1 

18  18-- 

6    0  6                  1 

181  l-~ 

<FR- -67                  1 

1  8  1  2  -  - 

3    18                  1 

1 8 13  — 

■no 23                1 

1  8  1  4 

+    '3  3                  1 

1815  — 

9    11                   1 

181  6-- 

3    O  3                  1 

1  d  1  { 

,;FR-...-67                  1 

1813  — 

9    11                   1 

1  8  1  9  -  - 

:,;T0 23                  1 

1  8  2  0  -  ~ 

4-    3  3                  1 

1  8  9  1  -  -- 

9   i  i                1 

1  C'  *"*  ■~|  — .  .„ 

1  L1  U.  L. 

4    04                  1 

1  8  2  3  -  - 

:to 44                1 

1824  — 

3/R---77                 1 

1  i-i  •-  cr 
J.  O  C  •_' 

1   01         ■   1 

1  8  d  6  — 

7    07                  1 

1  O  il  i' 

C  I.. !::! 37                  1 

1  '-'  '"'  o  ...  _ 

1  '—'  u—  O 

UP 27'                 1 

1  8  2  9  — 

XFR 67                  1 

1  8  3  0  -  - 

2    0  2                  1 

1831  — 

5    0  5                  1 

1  O  O  j 

X  =  Y 5  0                  1 

1 8  3  3  -  - 

1   01           1 

1834  — 

8 10                 1 

1  i_i  "I  c 

1  o  •  ".•  •_'  —  — 

6    0  6                  1 

8  3  6  -  -    8    10                   1 

I  8  90--  X  F  R •■  6  7 

o  o  "7 '.•■  ir  p c  -? 

O  C<   f     ft  1  K      0  i 

.891--    7    07 

8  38--  7  07        : 

.892--    4    04 

3  39--    7    07                   : 

.  l-'  -•  ■ J  """           —  —  2%  i-f 

848  — GTG 44 

894--  Y  TO 4  8 

t~i  ■'  i                  i~*   .*  f**i          "*?  ""•                   • 

L 895 —    7  ' 07 

842--     1 01 

L  8  9  6  -  -    6    0  6 

8  4  3  -  -    0    0  0 

o  q  7 •■.■•  r  p c  "7 

844  — XFR 67 

.  8  9  8  -  -  6  0  6 

845--    7    07 

'899--  2 82 

846--    3    0  3 

1 9  0  0  —  H I  V 3  5 

8  4  7  —  G  T  0 4  4 

[901—    G    15 

848  — S/R 77 

I  9  O  2  -  -  Y  TO 4  0 

849--    1    01 

i  9  0  3  -  -    8 0  8 

850--    0  00 

I  9  0  4  -  -  C  L  X -  3  7 

851—    DM 25 

L905  — X=Y 50 

p  cr  •;. _j_   .._ ■-■  ■-. 

L  9  O  6  —    1    - O  1 

853- -XFR 6  7 

1  9  O  7  -  -    9    11 

354--    7    07 

i  9  8  8  -  -    3    1  0 

855--    8    1  0 

L  9  0  9  -  -    3    18 

856--XFR 67 

1 9 1  0  -  -    9    11 

857--    X    3  6 

191 1--    5    -85 

3  58--    7    0  7 

1  9  1  2  -  -  •/,  T  0 2  3 

859~-    7    0  7 

1  9  1  3  -  -    o. 1  3 

8  60--    +    3  3 

191 4--    8    -- 10 

861— Y  TO -40 

1  9  1  5  -  -    9    -  1  1 

8  62--    7    07 

1 9 1 6--XT0-1 23 

8  6  3  -  -    5    85 

1  9  1  7  -  -    b     - 1  4 

364  — G  TO 4  4 

1 9 1 8  —  X  F  R -  6  7 

8  b  5  -  -  b  ••••  R ?  i 

1  9  1  9  -  -    8    0  8 

3  66--     1    0  1 

1920—    UP 27 

367--    7    ;-  0  7 

1921--XFR 6  7 

8  6  3  -  ~  >',  F  F: 6  7 

1  9  2  2  -  -    9    1  1 

O  t'  r1       {                  —  0  1 

1  9  2  3  -  -    0    80 

8  7  0  -  -    4    -  O  4 

1  ?  c.  4  —  ,•••,  /  1 J  O 

871- -CHS 32 

1925--    1     01 

872--    UP 27 

J  926--    9    11 

873--XFR 67 

1  9  2  7  -  -    3    ~  - 8  3 

874--    7    07 

1  9  2  8  -  -    8 18 

O  "7  6^ "7   „  |7j  "7 

1 9  2  9  -  -  E  E  X 2  6 

876  — FMT  — -42 

1  9  3  0  -  -    2    0  2 

877--     1     01 

1  9  3  1  -  -  X  T  0 -  2  3 

373--    UP 2  7 

1932--    o.    13 

879--FMT 42 

1  9  3  3  -  -  G  T  0 4  4 

8  80--    1    O  1 

1934--     1     ---01 

881—    4    04 

1935--    9    11 

-i  -i  •-!       cr        i-j  cr 
O  Oil  —   •_' U  •_' 

1  9  3  6  -  -    7    6  7 

883--    8 00 

1937--    1     81 

8  8  4  ~  -  X  T  0  -  - 2  3 

1938--XFR 67 

p  ri  cr .   ._.  _  _  ■)  •-. 

1  939-- I  NH 31 

8  8  6  -  -  X  F  R 6  7 

194  0--.    t,     14 

8  8  7  -  -    6    O  6 

1941 --    UP 27 

888--     I    01 

1942--XFR 6  7 

889--    UP 2  7 

1  9  4  3  -  -    8    80 
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1944  — X>Y 53 

1998  — FMT  — 

-42 

1945--  1 01 

1 999--  1  - - 

-01 

1946  —  9  11 

2006--  UP- - 

Cm     1 

1947  —  7 07 

2001— FMT  — 

-42 

1943--  1  01 

2  802--  1  -- 

-81 

1  9  4  9  -  ~  i  01 

2  0  0  3  -  -  4 

-84 

1  9  5  9  -  -  X  T  0 2  3 

2  8  8  4-- X F  R - - 

-67 

1951-- 34 

2  8  8  5  -  -  7  - - 

-87 

1952--  b  14 

2  6  0  6  -  -  7 

-87 

1  9  5  3  —  5  85 

2007 - -  U P  — 

—  ti  f 

1954  — XTQ 23 

28  0  8  - - X  F  R - - 

-67 

1  9  5  5  -  -  -  3  4 

2  0  0  9  -  -  6 

-  8  6 

1  9  5  6  -  -  o.  1  3 

20 10--  3  — 

-83 

195?--  b  14 

201  1--X>Y  — 

—  c  •"' 

1958—  UP 27 

28 12--  1  — 

-01 

1  9  5  9  ~  -  8  1  0 

28 13--  3  -- 

-8  3 

1  9  6  0  -  -  1  01 

28 14--  4  -- 

-  8  4 

1 9  6 1  -  -  X  =  Y 58 

•-i  tj  j  r   _   cr   

-85 

1  9  6  2  -  -  1  81 

2016— END  — 

-4  6 

1 9  6  3  -  -  3  1  1 

1  9  6  4  -  -  7  87 

1  9  6  5  ~  -  1  - -  0  1 

1  9  6  6  -  -  G  T  0  ~ 4  4 

1  9  6  7  -  -  1  0] 

1968  —  9  -- 11 

* 

1  9  6  9  -  -  3  8  3 

1  9  7  0  -  -  8  1  0 

19  71  —  8  0  0 

1972---  UP 2  7 

1 9  7  3  —  X  P  R 6  7 

1  9  7  4  -  -  7  0  7 

1  9  7  5  -  -  6  0  6 

1  9  7  6  -  -  X  >  Y 5  3 

1  9  7  7  -  -  1  01 

1978--  9 11 

1  9  7  9  -  -  8  1  0 

1  9  8  0  -  -  8  1  0 

1981--EEX- 26 

1  9  8  2  -  -  2  8  2 

1 9 8 3  - -  UP 27 

193  4--.  o.  1  3 

I  9  g  c;  _  _  _ •:,  ^ 

1  9  8  6  -  -  Y  T  0 4  0 

I  g  g  7  —  ,-,  j  o 

1988- -CHS  — -32 

■ 

1  9  8  9  -  -  2  8  2 

1  9  9  0  -  -  4 04 

1  9  9  1  -  ~  U  88 

1992  —  .UP 27 

1  9  9  3  -  -  a 1  3 

1  q  9  4  -  -  + 33 

1995--XFR 6  7 

1996  —  7  0  7 

1  9  9  7  -  -  7  8  7 
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c 

C      PEAPS  (FORTRAN  VERSION) 
C 

DIMENSION   HEADC20),  0(15),  T(15),  SALINU5),  S(15) 

DIMENSICN   SGRAD(15),  C0(6),  C S ( 6  ) ,  PRC13),  X(8) 

DIMENSION   IH(6) 

INTEGER  *  2   IBUF(61,  121) 

DATA  X/.  12,. 2 5, .3 8,. 52, .6  7, .84, 1 .06, 1.28/ 

INTEGER  *  2   IB/'   '/ 

INTEGER  *  2   IQOT/«.  »/ 

INTEGER  *  2   ISTAR/**  «/ 

PI  =  3.  141592 

RAD  =  PI  /  180.0 

90  DO  91  I  =  1,  6 
CD( I)  =  0.0 

91  CS(  I)  =  0.0 
C 

C      READ  AND  PRINT  HEADING  CARD  (80  COLUMNS  MAX) 
C 

READ  (5,  900,  END  =  9999)  HEAD 

900  FORMAT  (20A4) 
WRITE  (6,  901)  HEAD 

901  FORMAT  (1H1,  20A4) 
C 

C      READ  AND  PRIMT  INPUT  DATA  CARDS 
C 

READ  (5,  902)  NPTS,  MODE,  IBT,  IWH,  ISS,  IF  CM,  SD,  RD, 
-f  SF 

902  FORMAT  (615,  3F10.2) 
FCM  =  FLOAT ( IFOM) 

SS  =  FLOAT! ISS) 

WRITE    (6,     903)     NPTSV     MODE,     I6T,     IWH,     ISS,     IFOM,     SD, 

*  RD,  SF 

903  FORMAT  (1H0,  22HNUMSER  OF  SSP  POINTS  =,  15// 

*  19H  RAY  PLOTTER  MODE  =  ,  15  //14H  BOTTOM  TYPE  =, 

*  15  //14H  WAVE  HEIGHT  ^,  15,  3H  FT  // 

*  15H  SYSTEM  SIGMA  =,  15  //6H  FOM  =,  1 5  // 

*  13H  SOURCE  DEPTH  ,  F10.2,  3H  FT  // 

*  15H  RECEIVER  DEPTH  ,  F10.2,  3H  FT  // 

*  17H  SOURCE  FREQUENCY  ,  F10.2,  3H  HZ  ) 
WRITE  (6,  901)  HEAD 

WRITE  (6,  904) 

904  FORMAT  C 1H0, 5X, 5HP0INT f IX t 5HDEPTH ,4Xf  1 1HTEMPERATURE , 

*  IX,  8HSALINITY,  IX,  11HS0UND  SPEED,  IX, 

*  8HGRADIENT) 
C 

C      SPEED  OF  SOUND  PROGRAM 
C 

DO  102  I  =  1,  NPTS 

READ  (5,  905)  D(I),  T(I),  SALIN(I) 

905  FORMAT  (3F10.2) 

WRITE  (t>,  906)  I.  D(I),  Til),  SALIN(I) 

906  FORMAT  (1H0,  110,  3F10.2) 

T(I)  =  ( T (  I  )  -  3  2.0)  *  5.0  /  9.0 

D(l)  =  D( I )  *  0.3C48 

02  =  D(  i)     *  D( I ) 

D4  =  D2  *  D2 

T  2  =  T  (  I )  *  T  (  I  ) 

T3  =  T2  *  T(I) 

51  =  SALIN( I )  -  3  5.0 

52  =  SI  *  SI 

S(I)  =  1449.3  t  1.5848E-2  *  0(1)  *  Ir572E-7  *  D2 

*  -  3.46E-16  *  04  t  4*537  *  T(I)  -  5.356E-2  *  T2 

*  +  2.604E-4  *  T3  ♦  1.19  *  SI  +  9.6E-2  *  S2 

*  +  1.53<tE-4  *  T2  *  D(I)  -  7.19E-9  *  T(I)  *  D2 

*  -  1.2E-2  *  SI  *  TC21 
S(I)  =  S( I )  /  0.3048 

D(I)  =  DC  I )  /  0.3048 
WRITE  ( 6,  907)  SCI) 

907  FORMAT  (1H-,  40X,  F10.2) 
IF  (I  .GT.  1 )  GO  TO  100 
SGRAD( 1 )  =  0.0 
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GO  TO  101 

100  J  =  I  -  1 

SGRAD(I)  =  (S(I)  -  S(J1J  /  (0(1)  -  DiJ)) 

101  WRITE  (6,  908)  SGRADC I ) 
908  FORMAT  (1H0,  50X,  F10.4) 

102  CONTINUE 
C 

C      SPEED  OF  SOUND  PROFILE  CRITICAL  POINTS  PROGRAM 
C 

I  =  0 

I  FLAG  =  0 

IRUN  =  -1 

J  =  1 

K  =  2 

NG  =  0 

NTEST  =  -1 

103  I  =  I  +  1 
J  =  J  +  I 

NTEST  =  NTEST  +  1 

104  I  RUM  =  IRUN  *  1 

IF  (J  .GT.  NPTS)  GO  TO  113 
IF  (S( I )  .EQ.  S( J) )  GO  TO  103 
IF  (NTEST  .EQ.  0)  GO  TO  105 
NTEST  =  0 

cs(K)  =  sen 

C  D ( K )  =  0(1) 

K  =  K  +  1 

NG  =  NG  f  1 

IF  (S( J  )  .LT.  S( I  J )  GO  TO  105 

I  FLAG  -  1 

105  IF  (IFLAG  .EQ.  1)  GO  TO  110 

106  IF  (S(J)  .LT.  Sili)  GO  TO  108 
IF  (IRUN  .EQ.  OJ  GO  TO  107 
CS(K)  =  S< I) 

CD(K)  =  Oil) 
K  =  K  +  1 
NG  =  NG  +  1 
IFLAG  =  1 
GO  TO  104 

107  IFLAG  =  1 
GO  TO  104 

108  IF  (IRUN  cGT.  0)  GO  TO  109 
CS(K)  =  SC I) 

CD(K)  =  D(  I) 
K  =  K  *  1 
NG  =  NG  +  1 

109  IRUN  =  IRUN  *  1 
1  =  1  *  1 

J  =  J  +  1 

IF  (J  .GT.  NPTS)  GO  TO  112 

GO  TO  106 

110  IFLAG  =  0 

IF  (S( J)  .GT .  S( I ) )  GO  TO  111 

C  S  (  K )  =  S  ( I ) 

CD(K)  =0(1) 

K  =  K  +  1 

NG  =  NG  «■  1 

GO  TO  104 

111  I  ■  I  +  1 
J  =  J  *■  1 

IRUN  =  IRUN  +  1, 

IF  (J  .I.E.  NPTS)  GO  TO  110 

112  CS(K)  =  S( I) 
CD(K)  =0(1) 
NG  =  NG  +  1 

113  CS<  1)  =  S(l) 
CD( 1)  =  0(1) 
CS( 6)  =  S(NPTS) 
CD (6)  =  D(NPTS) 

C 

C      SCUND  SPEED  PROFILE  PROGRAM 

C 
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WRITE  (6.  901)  HEAD 
WRITE  (6,  909} 

909  FORMAT  OHO,  35HSOUND  SPEED  PROFILE  CRITICAL  POINTS  // 
*  5Xi  6H  DEPTH,  4Xf  14HSPEED  OF  SOUND  ) 

DO  114  1  =  It  6 

114  WRITE  (6,  910)  CD(I),  CS(I) 

910  FORMAT  (IHOt  2F10.2) 
DO  115  I  =  It  61 

DO  115  J  =  1,  121 

115  IBUF( I,  J)  =  IB 
WRITE  (6,  901)  HEAD 
WRITE  (6,  911) 

911  FORMAT  (IHOt  37HSPEED  OF  SOUND  (FT/SEC)  VS  DEPTH  (FT)) 
SMIN  =  CS(1) 

SMAX  =  CS(6) 

DO  116  I  =  1,  6 

IF  (CS( I)  .EQ.  0.0)  GO  TO  116 

IF  (CSU)  .LT.  SMIN)  SMIN  =  CS  (  I  ) 

IF  (CS(I)  *GT.  SMAX)  SMAX  =  CS ( I ) 

116  CONTINUE 

PR  I  =  (  SMAX  -  SMIN)  /  12.0 
PR«1)  =  SMIN 

V  s  PR (  1 ) 

DO  117  I  =  2,  13 

V  =  V  +  PRI 

117  PR( I)  =  V 

WRITE  (6.  912)  (PR(I),  I  =  1,  13) 

912  FORMAT  (F8.1,  12F10.1) 
DO  121  I  =  It  61 

IF  (MOD  KI-1),  10)  .EQ.  0)  GO  TO  119 
DO  118  J  =  It  121 »  10 

118  IBUFt I,  J)  =  IDOT 
GO  TO  121 

119  DO  120  J  =  1,  121 

120  1BUF( I,  J)  =  IDOT 

121  CONTINUE 

WRITE  (6,  913)  (IBUF(3»  J)t  J  =  It  121) 

913  FORMAT  (1H  t  5X,  121A1) 

PRI  =  i PR (13)  -  PRl  1) )  /  120.0 

PPDI  =  D(NPTS)  /  60.0 

J  =  NPTS  -  1 

DO  123  I  =  1,  J 

PPD  =  0(1) 

PPR  =  S(I ) 

PPRI  =  SGRAD(I  +  1)  *  PPDI 

122  IPPD  =  (PPD  /  PPDI)  *  1.5 
1PPR  =  (PPR  -  PR( 1) )  /  PRI 
IBUF(IPPD,  IPPR)  =  ISTAR 
PPD  =  PPD  +  PPDI 

IF  (PPD  .GT.  D( I  +  1))  GO  TO  123 
PPR  =  PPR  +  PPRI 
GO  TO  122 

123  CONTINUE 

MND   =    D(NPTS)    /    6.0 

MNDP    =    0 

DO    126     I    •-=    1,    61 

IF    (I    .  EQ.    1)    GO    TO    125 

IF    (MOD(U-l),     10)     .EQ.    0)     GO    TO    124 

WRITE    (6,     913)     (IBUFdi     J),     J    =     It     121) 

GO    TO    126 

124  MNDP  =  MNDP  +  MND 

125  WRITE  (6,  914)  MNDP,  (IBUF(I,JJ,  J  =  1,  121) 

914  FORMAT  (1H  ,  15,  121A1) 

126  CONTINUE 
C 

C      RAY  TRACE  PROGRAM 


C 


I  FLAG  -    0 

IF  (MODE  .NE.  2)  GO  TO  12  7 


C 

C      READ  IN  MODE  2  RAY  TRACE  EMISSION  ANGLES 

C 
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READ  (5,  905)  THMIN,  QTH,  THMAX 

127  IF  (MODE  .EQ.  1)  GO  TO  134 
DO  128  I  =  1»  61 

DO  128  J  =  1,  121 

128  IBUF  ( I ,  J)  =  IB 
WRITE  (6,  901)  HEAD 
WRITE  (6,  915) 

915  FORMAT  (1H0,  24HRANGE  (NM)  VS  DEPTH  (FT)  ) 
PRI  =  5.0 

V  =  -5.0 

DO  129  I  =  1,  13 

V  =  V  *■  PRI 

129  Pk(  I)  =  V 

WRITE  (6,  912)  (PR(I)t  I  =  1,  13) 

DO  133  I  =  It  61 

IF  (MOD  ((I  -  1),  10)  .EQ.  0)  GO  TO   131 

DO  130  J  =  1»  121,  10 

130  IBUF( I,  Ji  =  IDOT 
GO  TO  133 

131  DO  132  J  -  1,  121 

132  IBUF( I,  J)  =  IDOT 

133  CONTINUE 

PRI  =  3000.0 

RMAX  =  360000.0 

WRITE  (6,  913)  (IBUF(3,  J),  J  =  1,  121) 
C 
C       CONVERGENCE  ZONE  TEST 


C 


IF  (MODE  .EQ.  2)  GO  TO  219 

134  IF  (CS(6)  .GT.  CS(2))  GG  TO  136 

135  WRITE  (6,  916) 

916  FORMAT  (1H0,  19HN0  CONVERGENCE  ZONE  ) 
ICZ  =  0 
GO  TO  139 

136  M  =  2 

137  M  =  M  +  1 

IF  (CS(M)  .EQ.  0.0)  GO  TO  135 
IF  (CS(2)  .GT.  CS(M)}  GO  TO  138 
IF  (M  .EQ.  6)  GO  TO  135 
GO  TO  137 

138  ICZ  =  1 

139  M  =  0 

200  M  =  M  +  1 

IF  (  SD  .GE.  D(M)  )  GO  TO  200 

SV  =  ( ( SD  -  D(M  -  1))  *  SGRAD(M))  +  S(M  -  1) 

THP  =  (2.0  -  (2.0  *  SV  /  CS(2)))  **  0.5 

THPP  =  (ABS  (2.0  -  (2.0  *  SV  /  CS(6))))  **    0.5 

IF  (ICZ  .EQ.  0)  GO  TO  227 

THMIN  =  THPP 

I  PLOT  =  1 

ISR  =  1 

201  IF  (MODE  .EQ.  1)  IPLOT  =  0 
TH  =  THMIN 

Z  =  SD 
R  =  0.0 

202  M  =  0 

20  3  M  =  M  +  1 

IF  (D(MJ  .GT.  Z)  GO  TO  204 
IF  (M  .tQ.  NPTS)  GO  TO  208 
GO  TO  203 

204  IF  (Z  .LT.  0.0)  GO  TO  208 

205  RADIUS  =  (-SGRAD(M)  *  (1.0  -  i ( TH  *  TH)  /  2.0)))  / 
*  (  (  SGR  AD  (  M  )  *  (  Z  -  D(  M  -  1 )  )  *  S  ( M  -  I }  )  ) 

DR  =  5.0  *  (<D(M)  -  D(M  -  1))  **  0.5) 

206  RDR  =  DR  *  RADIUS 
r  —  r  +•  n n 

Z  -    I    +  MRADIUS  *  (DR  *  DR)  )  /  2.0)  *  (  TH  *  DR ) 

TH  =  TH  +  RDR 

IF  ( IPLOT  .EQ.  0)  GO  TO  207 

IPPD  -  (Z  /  PPD1  )  «■  1.5 

IF  (IPPD  .GT.  61 )  GO  TO  207 

1 FPR  =  R  /  PRI 


157 


IF  ( IPPR  .GT.  121  )  GO  TO  207 
IBUFUPPD,  IPPR)  =  I  STAR 

207  IF  (Z  .LT,  D(M  -  1))  GO  TO  202 
IF  (D(M)  .LT.  I)    GO  TO  202 

IF  (RMAX  .GE.  k)     GO  TO  206 

208  GO  TO  (209,  212,  213,  216,  217,  222,  220) 

209  IF  (D(NPTS)  .LT.  Z)  GO  TO  210 
IF  (IFLAG  .EQ.  1 )  GO  TO  211 
THMIN  =  THMIN  +  (0.5  *  RAD) 

I  PLOT  =  1 
ISR  =  1 
GO  TO  201 

210  THMIN  =  THMIN  -  (0.1  *  RAO) 
I  PLOT  =  0 

IFLAG  =  1 
ISR  =  1 
GO  TO  201 

211  THBB  =  THMIN 
IFLAG  =  0 
ICZ  -  1 

I  PLOT  =  1 
CZMAX  =  R 
ISR  ~    2 
GO  TO  201 

212  THMIN  =  THP 
I  PLOT  =  1 
ISR  =  3 

GO  TO  201 

213  IF  (R  .GT.  (CZMAX  /  2.0))  GO  TO  214 
THMIN  =  THMIN  *■  (0.1  *  RAD) 

IPLOT  =  0 
ISR  =  3 
GO  TO  201 

214  IF  (R  .GT.  CZMAX)  CZMAX  =  R 

215  CZMIN  =  R 

THMIN  =  THMIN  <■    (0.5  *  RAD) 
IPLOT  =  1 
ISR  =  4 
GO  TO  201 

216  IF  (R  .LT.  CZMIN)  GO  TO  215 
THMIN  =  THMIN  -  (0.4  *  RAD) 
IPLOT  =  0 

ISR  =  5 
GO  TO  201 

217  IF  (R  .GT.  CZMIN)  GO  TO  218 
CZMIN  =  R 

THMIN  =  THMIN  *-  (0.1  *  RAD) 
ISR  =  5 
GO  TO  201 

218  THMIN  =  THMIN  -  (0.1  *  RAD) 
IPLOT  =  1 

ISR  =  6 
GO  TO  201 

219  IF  (IFLAG  .EQ.  1)  GO  TO  222 
ISR  =  7 

GO  TO  201 

220  THMIN  -~  THMIN  +  DTH 

IF  (THMIN  .LE.  THMAX)  GO  TO  219 
ISR  =  1 
GO  TO  223 

221  MODE  =  1 
GO  TO  134 

222  IBB  =  1 

CZMIN  =  CZMIN  /  3.0 
CZMAX  =  CZMAX  /  3.0 
I  PLOT  =  1 
ISR  =  2 
2  23  MNDP  =  0 

DO  22  6  I  =  1,  61 

IF  ( I  .EQ.  1)  GO  TO  225 

IF  (MODHI-i),  10)  .EQ.  0)  GO  TO  224 

WRITE  (6,  913)  (IBUFCIf  J),  J  =  1,  121/ 
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24 

25 
26 


C 
C 

c 


227 


17 


28 


229 


30 
31 


33 


GO  TO  226 

MNDP  =  MNDP  +  MND 

WRITE  (6,  914)  MNDP, 

CONTINUE 

GO  TO  (  221,  227) ,  ISR 

CONTINUE 


f  I BU  F  (  I  •  J ) 


1  ,  121) 


TRANSMISSION  LOSS  PROGRAM 


B  =  FOM 
IF(80.G 
IF(100. 
B  =  (B 
B  =  10. 
IR  =  B/ 
RMAX  = 
WRITE  ( 
WRITE  ( 
FORMAT 
DO  228 
DO  228 
IBUF( I, 
PRI  =  R 
PR(1)  = 

V  =  0,0 
DO  229 

V  =  V  4- 
PRU)  = 
IPFOM  = 
DO  233 
IF  (MOD 
DO  230 
IBUF( I, 
IF  (I  . 
CONTIUU 
GO  TO  2 
DO  232 
IBUF( I, 
IF  (1  . 
CCNTINU 
CCNTINU 
PRI  =  R 
PPDI  = 
IF (CSC  2 
IFCSD.G 
I  F  {  RD  .  G 


E.I  FOM)  B  =  80.0 

LE.IFOM)  B  =  100,0 

+  15.0)  /  20.0 

0**B 

20000.0 

(IR  *  1)  *  2  00  00 

6,  901)  HEAD 

6,  917) 

(lHOi  10HRANGE  (NMJ 

I  =  1,  61 

J  =  1,  121 

J)  =  IB 
MAX  /  24000.0 

0.0 


IX,  5HTLGSS  ,  4X,  4HPDET) 


I  =  2,  13 
PRI 
V 

( ( IFGM  -  40)  / 
I  =  1,  61 
UI-1),  10)  „EQ, 
J  =  1,  121,  10 
=  IDOT 
IPFGM)  I  BUF (I 


?) 
0) 


J) 

EQ 

E 

33 

J  =  1,  121 

=  IDOT 
IPFOM) 


GO  TO  231 
J)  =  I  STAR 


ISO 
IRAP 
RS  = 


J) 

EQ, 

E 

E 

MAX 

2.0 

)  -LT.CS( 1 )  )  GO 

T.CD(2>)  GO  TO 

T.CD(2) )  GO  TO 


I  BUM  If  J)  =  I  STAR 


/  120.0 


TO  304 

30? 
301 


=  1 
=  0 


301 


AR 

*  )/ 
RT  =  CO 
SSL  =  ( 
G  =  ABS 
DFL  =  7 

*  ( 
GO  TO  3 
IRAP  = 
RL2  =  C 
IFiSD.G 


02 

03 

04 
05 

06 
0  7 


308 


RL2  = 
GO  TO 
RL2  = 
GO  TO 
IF(RD.G 
GO  TO  3 
IRAP  = 
I  SD  =  0 
GO  TO  3 
IF(NC.F 
IF(  ICZ. 
GO  TU  3 
CZSTRT 


COS ( CS ( 1 ) /CS ( 2 ) ) * ( 2 . 0/3 . 0 ) * (CS ( i } *{ CD( 2 ) ~CD( 1 ) ) 

(CS(2i  -  CS(  1J ) 

(2)/(6.0*ARC0S(SV/CS<2) ) ) 

0 . 8/ RS ) - (  ( WH* S  F/ 1 000 . 0 ) *  *0 . 7 1 ) 

((CS(3)-CS(2) )/(CD(3)- -CD( 2)  )  ) 

44.0*<  (SF/10GO.O)**{-5.  0/3.0)  )*( G**( -1 .0/3.0)  )* 

CD (2)  **  (-3.0)) 

06 

1 

S( 2 )*(CD< 3 )-CDi 2) )/<CS(2)-CS(3J) 

7*CD( 2) )  GO  TO  30  2 

L2*<RD-CD(2))/4.5 

05 

L2*(SD-CD(2) J/4.5 

05 

T.C0(2) )  GO  TO  304 

01 

0 

07 

Q.l)  DFL  =  0.0 

EQ.l)  GO  TO  308 

10 

=  CZMIN 


CZSTOP  =  CZMAX 
IFURAP  .EQ.l)  GO  TO  310 
IF( ISD.EQ.l )    GO  TO  310 
ICSC  =  1 
M  =  1 
3  09  M  =  M  +  1 

IF(CS(M+1) ,LT.CS(MJ )  GO  TO  309 

RPT  =  (CZMIN/  8.0)*((CD(M)-CD12)  )/(SD-CD(2)  ))**0.5 

GO  TO  313 

310  FTEST  =  3  2.0  *  <(SV  /  SF)  **  0.7) 
IF(FTEST.GT.CD(2M  GO  TO  311 

GO  TO  312 

311  ISO  =  0 

312  ICSC  =  0 

313  IFC6.LT. IflT)  GO  TO  314 
IF(5.LT.I3T)  GO  TO  315 
IF<4.LT.IBT)  GO  TO  316 
IFO.LT.IBT)  GO  TO  317 
IF(2.LT.IBT)  GO  TO  318 
BT  =  0.66 

GO  TO  319 

314  BT  =  2.0 
GO  TO  319 

315  BT  =  1.41 
GO  TO  319 

316  6T  =  1. 18 
GO  TO  319 

317  BT  =  1.0 
GO  TO  319 

318  BT  =  0.85 

319  3LMA  =  30,0 
IF(3500.0.LT.SF }  GO  TO  320 
IF(2000.0.LT.SF)  GO  TO  321 
BLMA  =  50.0 

IFUOOO.O.LT.SF)  GO  TO  322 
IF(500.0.l.T.SF)  GO  TO  323 
IF(IOO.O.LT.SF)  GO  TO  324 
Bl  =  0.2 

BO  =  -0.83 
GO  TO  325 

320  Bl  -  0.35 
BO  =  7.12 
GO  TO  325 

321  Bl  =  0.28 
BO  =  8.2) 
GO  TO  325 

322  Bl  =  0.19 
BO  =  2.09 
GO  TO  325 

323  Bl  =  0. 17 
BO  =  1.97 
GO  TO  325 

324  Bl  =  0. 18 
BO  *  2*06 

325  SFKHZ  =  SF/1000.0 
SFKHZ2  =  SFKHZ*SFKHZ 

A  =  (  (0.1*SFKHZ2)/(  1.0+  SFKHZ2)  )  4-  ( ( 40* SFKHZ2 ) / ( 4 100.0 
*     +  SFKHZ2J ) 
A  =  A  /  1000.0- 
ICZRC  =  1 
H  =  1.0 

DRWRK  =  2000.0 
ORPERM  =  2000.0 
IF  (ICZ  .EQ.  1)  GO  TO  3251 
IF  (CS(6)  ,GT.  SV)  GO  TO  3250 
THP  =  ARCOSISV  /  CS(6) ) 

3250  THB3  =  THP 

3251  R  =  0.0 

THBL  =  PI  /  2.0 

326  IF(THBB.GT.ThBL)  H  .=  H  *  1 . 0 
R  =  R  *  DRWRK 

THBL  =  ATANf  (R*UCD(6  J-SDJ/3.G)  )/ 


1G0 


327 
328 

329 


330 


331 


332 

3 

333 
334 


<H*( ( 2.0*CD(6 )-SD)/3.0) ) ) 
TH6LA    =    THBL    *     ( 180*0    /PI) 
IF     (BLMA    .LT.     THBLA)     THBLA    =    BLMA 
BL    =    H    *    BT    *     ((THBLA    *    81)     «•    BO) 
STL    =     <20.0*ALGG10(R/C0S<THBL) )     + 

(  A*{R/CCS(THBD)  )     *    6L    -    6.0) 
IF( ISD. EQ.O)    GO    TO    330 
1SR    =    1 

A    «■    DFL)     *    R 
.LT.    R)     GO    TO    327 

L    +    20.0*ALOG10(R) 

^28 

L    +    10.0*ALOG10(R) 
.GT.    R)    GO    TO    329 

L    +    (SSL    *    R) 

BS(  STL    -    TLM    /     (-10.0) 

0.0    *    (ALCC-lOd.O    ♦    (10.0    **    F))) 

.LT.TL)     TL    =    STL 

TL    -    AO 

(330$     334,     :>;>o, 

P  .EQ.O)    GO    TO    334 

2 
=     (RL2    **    0,5)     *    (RS    /    2.0) 


TL 
IF 
TL 
GO 
TL 
IF 
TL 


=     ( 

(RT 
=    T 
TO 
=    T 
(RS 
=    T 
=     (A 
AD    =    1 
IF(STL 
STL    = 
GO    TO 
IFC  IRA 
ISR    = 
RAPTES 
I  F  ( R  .  L 
I  RAP    = 
GO    TO 
IF    (SO 


+    10.0*ALOG10(RT) 


338,     339),     ISR 


=    2 

R 
TO 
=    2 

R 

TO 


335 


336 


337 
338 


339 


TL 

GO 
TL 

GO 
IF(  ICZ 
ISR    = 
TL    =    2 
IFU.L 
IF(R.G 
1 FC  ICZ 
DRWRK 
R    =    CZ 
ICZRC 
GO    TO 
RANGE 
GCZMAX 
GB    -    0 
RP    =    R 
GCZ    = 
TL    =    T 
GO    TO 
CZSTRT 
CZSTOP 
DRWRK 
ICZRC 
GO    TO 
IF ( IDS 
ISR    = 
TL    *    1 
GO    TG 
PPD    = 
PPR    = 
IPPD    = 
IF    (IP 
I  PPR    = 
IF     (IP 
IF    (  IP 
IBUF(  I 


T. RAPTES)     GO    TO    331 

0 
3  30 

.LT.RDJ    GO    TO    332 
0.0*ALOG10( (SD*SD    t 

*    R)    **    0.5) 


R*R)**0.5) 
R*R)**0.5) 


*  A*((SD*SD    + 

*  A*((RD*RD    + 


0.0*AL0G10( ( RD*RD    f 

*    R)    **    0.5) 
329 

.EQ.O)     GO    TO    338 
3 

0.0*AL0G10(R)     *    A*R 
T. CZSTRT)     GO    TO    337 
T. CZSTOP)     GO    TO    336 
RC.E0.2)     GO    TO    335 
=     (CZSTOP    -    CZSTRTI/10.0 
STRT    -    DRWRK 
=    2 
3  26 
=    R    -    CZSTRT 

=    15.5    -    0.0005*<CZSTGP    -    CZSTRT) 
.36    /    (GCZMAX*GCZMAX) 

ANGE*(GB**.5)    /     I  0. 1 8*( CZSTOP    -    CZSTRT)) 
(RP/GB)     *    EXPC-tRP*RP)    /     (2.0*GB)i 
L    -    GCZ 
337 

=    CZSTRT 

=  CZSTOP 
=  OR PERM 
=  1. 
329 
C.EQ.O)  GO  TO  339 

6.0*AL0G10(Ri  ♦  10.0*AL0G10(RPT)  +  A*R 

329 

STL 

R 

( ( PPD  -  40.0)  /  PPDI )  *  1.5 
PD  .GT«  61 )  IPPD  =  61 

(PPR  /  PRI )  +  1.0 
PR  .GT.  121 )  I PPR  =  121 
PD  .GT.  60)  IPPD  s  60 
PPD,  IPPR)  =  1STAR 


<•  CZfUN 
i-    CZMAX 


C 


PROBABILITY  OF  DETECTION  PROGRAM 


PD  =  0.5 
SE  =  FOM 


STL 


61 


340 


341 
34  2 

343 


918 


919 


344 


345 
346 

347 

9999 


XX  ■■= 

IF(XX 
PD  = 
I  =  8 
I  F  (  1 « 

IFiXX 
PD  = 
I  =  I 
IF  (  I . 

GO  TO 
PD  = 
I  F  (  SE 
PD  = 
PPD  = 
IPPD 
IF  (I 
IBUFC 
PP  = 
WRITE 
FORMA 
I F  ( RM 
WRITE 
WRITE 
FORMA 

WRITE 
WRITE 
MNDP 

MND  = 


DO 
IF 
IF 
IF 


34 
(  I 
(M 
(  I 


WRITE 

GO  TO 

IF  (I 

MNDP 

MND  = 

MNDP 

WRITE 

CGNTI 

GO  TO 

STOP 

END 


ABS(SE/SS) 
.EO.0.0)  GO  TO  343 
0,95 

64.LE.XX)  GO  TO  341 
,GT.X( I ) )  GO  TC  342 
PD  -  0.05 

-  1 
EQ.OJ  GO  TO  342 

340 
1.0 

.GT.O.O)  GO  TO  343 
1.0  -  PD 

61.0  -  (PD  *  10.0) 
=  PPD  +  0.5 

PPD  .GT.  61)  IPPD  =  61 
IPPD,  IPPR)  s  ISTAR 
R  /  2000.0 

(6,  918)  RP,  STL,  PD 
T  (1H0,  3F10.2) 
AX.GT.R)  GO  TO  326 

(6,  901)  HEAD 

(6,  919) 
T  (1H0,  35HTL0SS  (DB)  ABOVE  PDET  (PERCENT)  VS  , 
10HRANGE  (MM)  ) 

(  6t  912)  (  PR(  I)  ,  I  =  1,  13) 

(6,  913)  (IBUF(3,  J),  J  =  i,  121) 
=  40 

20 
7  I  =  1,  61 

.EQ.  1)  GO  TO  346 
0D((I-1),  10)  .EQ.  0)  GO  TO  345 

.GT.  50)  GO  TO  344 

(6,  913)  ( IBUF( I,  J) ,  J  =  1,  121) 

347 

.GT.  52)  GO  TO  345 
=  100 

-10 
=  MNDP  +  MND 

(6,  914)  MNDP,  (IBUF(I-tJ)f  J  =  It  121) 
NUE 

90 
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APPENDIX  E 
PEAPS  STORAGE  LOCATIONS  IN  THE  HP 9 810 A  CALCULATOR 

A.  The  operator  enters  all  angles  in  degrees.   The  zero 
degrees  direction  for  all  angles  is  horizontal  to  the  right. 
Clockwise  angles  (down)  from  the  horizontal  are  to  be  entered 
as  positive  (+)  angles.   Counter-clockwise  angles  (up)  from 
the  horizontal  are  to  be  entered  as  negative  (-)  angles. 

B.  All  angles  are  stored  and  used  in  radians.   Entered  angles 
are  automatically  converted  from  degrees  to  radians  for 
storage  by  the  program: 

degrees  =  -^-^r   *  radians 

loU 


a  4  180  *  j 

radians  =  *  degree: 


C.   Enter  all  temperatures  in  degrees  fahrenheit  (°F) 
Temperatures  are  automatically  converted  to  degrees 
centigrade  (°C)  for  use  in  the  SSP  program: 


°C  =  I  (°F  -  32) 


°F  =  32  +  (I-  *  °C) 


16  3 


D.  Enter  depths  in  feet  (ft) .   Depths  are  stored  in  feet 
and  are  converted  to  the  required  units  as  needed  in  the 
programs: 

meters  =0.304  8  *  feet 

feet  =  meters/0.304  8 

yards  =  feet/3.0 

nautical  miles  (NM)  =  yards/2000  =  meters/1852. 

E.  Ranges  are  in  feet  in  the  ray  trace  program  and  are 
automatically  converted  to  units  of  yards  for  use  in  the 
transmission  loss/probability  of  detection  program. 

F.  The  following  pages  of  program  forms  indicate  the 
parameter  storage  locations  used  in  the  PEAPS  model.   The 
terms  listed  are  defined  in  the  Table  of  Symbols  and 
throughout  the  body  of  the  text. 
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